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The physics program in Hall A at Jefferson Lab commenced in the summer of f 997 with a de- 
tailed investigation of the 16 0(e, e'p) reaction in quasielastic, constant (q,co) kinematics at Q 2 ~ 
0.8 (GeV/c) 2 , q ~ 1 GeV/c, and u) ~ 445 MeV. Use of a self-calibrating, self-normalizing, thin-film 
waterfall target enabled a systematically rigorous measurement. Five- fold differential cross-section 
data for the removal of protons from the lp-shell have been obtained for < p m i ss < 350 MeV/c. 
Six-fold differential cross-section data for < -Emiss < 120 MeV were obtained for < p m i ss < 340 
MeV/c. These results have been used to extract the Alt asymmetry and the Rl, Rt, Rlt, and 
Rl+tt effective response functions over a large range of i? m i ss and p m iss- Detailed comparisons 
of the lp-shell data with Relativistic Distorted- Wave Impulse Approximation (rdwia), Relativistic 
Optical-Model Eikonal Approximation (romea), and Relativistic Multiple-Scattering Glauber Ap- 
proximation (rmsga) calculations indicate that two-body currents stemming from Meson- Exchange 
Currents (MEC) and Isobar Currents (IC) are not needed to explain the data at this Q 2 . Further, 
dynamical relativistic effects are strongly indicated by the observed structure in Alt at p m iss ~ 300 
MeV/c. For 25 < -E m i ss < 50 MeV and p m i ss ~ 50 MeV/c, proton knockout from the lsi/2-state 
dominates, and romea calculations do an excellent job of explaining the data. However, as p m iss 
increases, the single-particle behavior of the reaction is increasingly hidden by more complicated 
processes, and for 280 < p m i ss < 340 MeV/c, ROMEA calculations together with two-body currents 
stemming from MEC and IC account for the shape and transverse nature of the data, but only about 
half the magnitude of the measured cross section. For 50 < -E m i ss < 120 MeV and 145 < p m i ss < 
340 MeV/c, (e, e'pN) calculations which include the contributions of central and tensor correlations 
(two-nucleon correlations) together with MEC and IC (two-nucleon currents) account for only about 
half of the measured cross section. The kinematic consistency of the lp-shell normalization factors 
extracted from these data with respect to all available ie O(e, e'p) data is also examined in detail. 
Finally, the Q 2 -dependence of the normalization factors is discussed. 

PACS numbers: 25.30.Fj, 24.70.+S, 27.20.+n 



I. INTRODUCTION 



Exclusive and semi-exclusive (e, e'p) in quasielastic 
(QE) kinematics [133] has long been used as a precision 
tool for the study of nuclear electromagnetic responses 
(see Refs. [1-4]). Cross-section data have provided in- 
formation used to study the singlc-nuclcon aspects of 
nuclear structure and the momentum distributions of 
protons bound inside the nucleus, as well as to search 
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for non-nucleonic degrees of freedom and to stringently 
test nuclear theories. Effective response-function separa- 
tions [134] have been used to extract detailed information 
about the different reaction mechanisms contributing to 
the cross section since they are selectively sensitive to 
different aspects of the nuclear current. 

Some of the first (e,e'p) energy- and momentum- 
distribution measurements were made by Amaldi et al. 
[5]. These results, and those which followed (see Refs. 
[1, 2, 6]), were interpreted within the framework of 
single-particle knockout from nuclear valence states, even 
though the measured cross-section data was as much as 
40% lower than predicted by the models of the time. The 
first relativistic calculations for (e, e'p) bound-state pro- 
ton knockout were performed by Picklesimer, Van Orden, 
and Wallace [7-9] . Such Relativistic Distorted- Wave Im- 
pulse Approximation (rdwia) calculations are generally 
expected to be more accurate at higher Q 2 , since QE 
(e, e'p) is expected to be dominated by single-particle in- 
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tcractions in this regime of four-momentum transfer. 

Other aspects of the structure as well as of the reaction 
mechanism have generally been studied at higher missing 
energy (E m i ss ). While it is experimentally convenient to 
perform measurements spanning the valence-state knock- 
out and higher E miss excitation regions simultaneously, 
there is as of yet no rigorous, coherent theoretical pic- 
ture that uniformly explains the data for all -E m iss and 
all missing momentum (p m i ss ). In the past, the theoret- 
ical tools used to describe the two energy regimes have 
been somewhat different. Miither and Dickhoff [10] sug- 
gest that the regions are related mainly by the transfer 
of strength from the valence states to higher E m i ss . 

The nucleus 16 O has long been a favorite of theorists, 
since it has a doubly closed shell whose structure is thus 
easier to model than other nuclei. It is also a convenient 
target for experimentalists. While the knockout of lp- 
shell protons from 16 O has been studied extensively in 
the past at lower Q 2 , few data were available at higher 
£" m i SS for any Q 2 in 1989, when this experiment was first 
conceived. 



A. lp-shell knockout 

The knockout of lp-shell protons in 16 0(e,e'p) was 
studied by Bernheim et al. [11] and Chinitz et al. [12] 
at Saclay, Spaltro et al. [13] and Leuschner et al. [14] 
at NIKHEF, and Blomqvist et al. [15] at Mainz at Q 2 
< 0.4 (GcV/c) 2 . In these experiments, cross-section data 
for the lowest-lying fragments of each shell were measured 
as a function of p m i ss , and normalization factors (relat- 
ing how much lower the measured cross-section data were 
than predicted) were extracted. These published nor- 
malization factors ranged between 0.5 and 0.7, but Kelly 
[2, 4] has since demonstrated that the Mainz data sug- 
gest a significantly smaller normalization factor (see also 
Table X). 

Several calculations exist (see Refs. [16-21]) which 
demonstrate the sensitivity [135] of the longitudinal- 
transverse interference response function Rlt and the 
corresponding left-right asymmetry Alt [136] to 'spinor 
distortion' (see Section IV A 1), especially for the removal 
of bound-state protons. Such calculations predict that 
proper inclusion of these dynamical rclativistic effects 
is needed to simultaneously reproduce the cross-section 
data, Alt, and Rlt- 

Fig. 1 shows the effective response Rlt as a function 
of Pmiss for the removal of protons from the lp-shell of 
16 O for the QE data obtained by Chinitz et al. at Q 2 — 
0.3 (GeV/c) 2 (open circles) and Spaltro et al. at Q 2 = 
0.2 (GeV/c) 2 (solid circles) together with modern rdwia 
calculations (see Sections IV and V for a complete dis- 
cussion of the calculations). The solid lines correspond to 
the 0.2 (GcV/c) 2 data, while the dashed lines correspond 
to the Q 2 = 0.3 (GeV/c) 2 data. Overall, agreement is 
good, and as anticipated, improves with increasing Q 2 . 
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FIG. 1: Longitudinal-transverse interference effective re- 
sponses Rlt as a function of p m i ss for the removal of protons 
from the lp-shell of ie O. The open and filled circles were ex- 
tracted from QE data obtained by Chinitz et al. at Q 2 = 0.3 
(GeV/c) 2 and Spaltro et al. at Q 2 = 0.2 (GeV/c) 2 , respec- 
tively. The dashed (Q 2 = 0.3 (GeV/c) 2 ) and solid (Q 2 = 0.2 
(GeV/c) 2 ) curves are modern rdwia calculations. Overall, 
agreement is good, and improves with increasing Q 2 . 



B. Higher missing energies 



Few data are available for 16 0(e, e'p) at higher -B m is S , 
and much of what is known about this excitation region is 
from studies of other nuclei such as 12 C. At MIT-Bates, a 
series of 12 C(e,e'p) experiments have been performed at 
missing energies above the two-nucleon emission thresh- 
old (see Refs. [22-26]). The resulting cross-section 
data were much larger than the predictions of single- 
particle knockout models [137]. In particular, Ulmcr et 
al. [23] identified a marked increase in the transverse- 
longitudinal difference St — Sl [138]. A similar increase 
has subsequently been observed by Lanen et al. for 6 Li 
[27], by van der Steenhoven et al. for 12 C [28], and most 
recently by Dutta et al. for 12 C [29], 56 Fe, and 197 Au 
[30]. The transverse increase exists over a large range 
of four-momentum transfers, though the excess at lower 
Pmiss seems to decrease with increasing Q 2 . Theoreti- 
cal attempts by Takaki [31], the Ghent Group [32], and 
Gil et al. [33] to explain the data at high -E m iss using 
two-body knockout models coupled to Final-State Inter- 
actions (FSI) have not succeeded. Even for QE kine- 
matics, this transverse increase which starts at the two- 
nucleon knockout threshold seems to be a strong signa- 
ture of multinuclcon currents. 
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FIG. 2: The experimental infrastructure in Hall A at Jeffer- 
son Lab at the time of this experiment. The electron beam 
passed through a beam-current monitor (BCM) and beam- 
position monitors (BPMs) before striking a waterfall target 
located in the scattering chamber. Scattered electrons were 
detected in the HRS e , while knocked-out protons were de- 
tected in the HRSh. Non-interacting electrons were dumped. 
The spectrometers could be rotated about the central pivot. 



II. EXPERIMENT 

This experiment [34, 35], first proposed by Bertozzi et 
al. in 1989, was the inaugural physics investigation per- 
formed in Hall A [36] (the High Resolution Spectrometer 
Hall) at the Thomas Jefferson National Accelerator Fa- 
cility (JLab) [37]. An overview of the apparatus in the 
Hall at the time of this measurement is shown in Fig. 2. 
For a thorough discussion of the experimental infrastruc- 
ture and its capabilities, the interested reader is directed 
to the paper by Alcorn et al. [38]. For the sake of com- 
pleteness, a subset of the aforementioned information is 
presented here. 



A. Electron beam 

Unpolarized 70 /iA continuous electron beams with en- 
ergies of 0.843, 1.643, and 2.442 GeV (corresponding to 
the virtual photon polarizations shown in Table I) were 
used for this experiment. Subsequent analysis of the data 
demonstrated that the actual beam energies were within 
0.3% of the nominal values [39]. The typical laboratory 
±4tT beam envelope at the target was 0.5 mm (horizon- 
tal) by 0.1 mm (vertical). Beam-current monitors [40] 
(calibrated using an Unser monitor [41]) were used to 
determine the total charge delivered to the target to an 
accuracy of 2% [42]. Beam-position monitors (BPMs) 
[43, 44] were used to ensure the location of the beam at 



TABLE I: The QE, constant (q, ui) kinematics employed in 
this measurement. At each beam energy, q ~ 1 GeV/c. 



-^bcam 


9e 


virtual photon 


0,,,, 


(GeV) 


n 


polarization 


(°) 


0.843 


100.76 


0.21 


0, 8, 16 


1.643 


37.17 


0.78 


0, ±8 


2.442 


23.36 


0.90 


0,±2.5, ±8, ±16, ±20 



the target was no more than 0.2 mm from the beamlinc 
axis, and that the instantaneous angle between the beam 
and the beamline axis was no larger than 0.15 mrad. 
The readout from the BCM and BPMs was continuously 
passed into the data stream [45]. Non- interacting elec- 
trons were dumped in a well-shielded, high-power beam 
dump [46] located roughly 30 m from the target. 



B. Target 

A waterfall target [47] positioned inside a scattering 
chamber located at the center of the Hall provided the 
H2O used for this study of 16 O. The target canister was 
a rectangular box 20 cm long x 15 cm wide x 10 cm 
high containing air at atmospheric pressure. The beam 
entrance and exit windows to this canister were respec- 
tively 50 ^m and 75 ^m gold-plated beryllium foils. In- 
side the canister, three thin, parallel, flowing water films 
served as targets. This three-film configuration was su- 
perior to a single film 3x thicker because it reduced the 
target-associated multiple scattering and energy loss for 
particles originating in the first two films and it allowed 
for the determination of the film in which the scattering 
vertex was located, thereby facilitating a better overall 
correction for energy loss. The films were defined by 2 
mm x 2 mm stainless-steel posts. Each film was sep- 
arated by 25 mm along the direction of the beam, and 
was rotated beam right such that the normal to the film 
surface made an angle of 30° with respect to the beam 
direction. This geometry ensured that particles originat- 
ing from any given film would not intersect any other film 
on their way into the spectrometers. 

The thickness of the films could be changed by vary- 
ing the speed of the water flow through the target loop 
via a pump. The average film thicknesses were fixed at 
(130 ± 2.5%) mg/cm 2 along the direction of the beam 
throughout the experiment, which provided a good trade- 
off between resolution and target thickness. The thick- 
ness of the central water film was determined by com- 
paring 16 0(e, e') cross-section data measured at q « 330 
MeV/c obtained from both the film and a (155 ± 1.5%) 
mg/cm 2 BeO target foil placed in a solid-target ladder 
mounted beneath the target canister. The thicknesses 
of the side films were determined by comparing the con- 
currently measured 1 H(e, e) cross section obtained from 
these side films to that obtained from the central film. In- 
stantaneous variations in the target-film thicknesses were 
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TABLE II: Selected results from the optics commissioning. 



parameter 


resolution 
(FWHM) 


reconstruction 
accuracy 


out-of-plane angle 


6.00 mrad 


±0.60 mrad 


in-plane angle 


2.30 mrad 


±0.23 mrad 


^/target 


2.00 mm 


±0.20 mm 


Ap/p 


2.5 x 10~ 4 





monitored throughout the entire experiment by continu- 
ously measuring the H(e, e) cross section. 

C. Spectrometers and detectors 

The base apparatus used in the experiment was a pair 
of optically identical 4 GeV/c superconducting High Res- 
olution Spectrometers (HRS) [48]. These spectrometers 
have a nominal 9% momentum bite and a FWHM mo- 
mentum resolution Ap/p of roughly 10~ 4 . The nom- 
inal laboratory angular acceptance is ±25 mrad (hor- 
izontal) by ±50 mrad (vertical). Scattered electrons 
were detected in the Electron Spectrometer (HRS e ), and 
knocked-out protons were detected in the Hadron Spec- 
trometer (HRSh) (sec Fig. 2). Before the experiment, 
the absolute momentum calibration of the spectrometers 
was determined to Ap/p = 1.5 x 1CP 3 [39]. Before and 
during the experiment, both the optical properties and 
acceptances of the spectrometers were studied [49] . Some 
optical parameters are presented in Table II. During the 
experiment, the locations of the spectrometers were sur- 
veyed to an accuracy of 0.3 mrad at every angular loca- 
tion [50]. The status of the magnets was continuously 
monitored and logged [45]. 

The detector packages were located in well-shielded de- 
tector huts built on decks located above each spectrome- 
ter (approximately 25 m from the target and 15 m above 
the floor of the Hall). The bulk of the instrumentation 
electronics was also located in these huts, and operated 
remotely from the Counting House. The HRS e detector 
package consisted of a pair of thin scintillator planes [51] 
used to create triggers, a Vertical Drift Chamber (VDC) 
package [52, 53] used for particle tracking, and a Gas 
Cerenkov counter [54] used to distinguish between ir~ 
and electron events. Identical elements, except for the 
Gas Cerenkov counter, were also present in the HRSh 
detector package. The status of the various detector sub- 
systems was continuously monitored and logged [45] . The 
individual operating efficiencies of each of these three de- 
vices was >99%. 



D. Electronics and data acquisition 

For a given spectrometer, a coincidence between sig- 
nals from the two trigger-scintillator planes indicated a 
'single-arm' event. Simultaneous HRS e and HRS/j singles 



events were recorded as 'coincidence' events. The basic 
trigger logic [55] allowed a prescaled fraction of single- 
arm events to be written to the data stream. Enough 
HRS e singles were taken for a 1% statistics 1 H(e, e) cross- 
section measurement at each kinematics. Each spectrom- 
eter had its own VME crate (for scalers) and FASTBUS 
crate (for ADCs and TDCs). The crates were managed 
by readout controllers (ROCs). In addition to overseeing 
the state of the run, a trigger supervisor (TS) generated 
the triggers which caused the ROCs to read out the crates 
on an event-by-event basis. The VME (scaler) crate was 
also read out every ten seconds. An event builder (EB) 
collected the resulting data shards into events. An ana- 
lyzer/data distributer (ANA/DD) analyzed and/or sent 
these events to the disk of the data-acquisition computer. 
The entire data-acquisition system was managed using 
the software toolkit CODA [56]. 

Typical scaler events were about 0.5 kb in length. Typ- 
ical single-arm events were also about 0.5 kb, while typi- 
cal coincidence events were about 1.0 kb. The acquisition 
deadtime was monitored by measuring the TS output-to- 
input ratio for each event type. The event rates were set 
by varying the prescale factors and the beam current such 
that the DAQ computer was busy at most only 20% of 
the time. This resulted in a relatively low event rate (a 
few kHz), at which the electronics deadtime was <1%. 
Online analyzers [57] were used to monitor the quality 
of the data as it was taken. Eventually, the data were 
transferred to magnetic tape. The ultimate data analy- 
sis was performed on the DEC-8400 CPU farm ABACUS 
[58] at the Massachusetts Institute of Technology using 
the analysis package espace [59] . 

III. ANALYSIS 

The interested reader is directed to the Ph.D. theses of 
Gao [60] and Liyanage [61] for a complete discussion of 
the data analysis. For the sake of completeness, a subset 
of the aforementioned information is presented here. 



A. Timing corrections and particle identification 

Identification of coincidence (e, e'p) events was in gen- 
eral a straightforward process. Software corrections 
were applied to remove timing variations induced by 
the trigger-scintillator circuit and thus sharpen all flight- 
time peaks. These included corrections to proton flight 
times due to variations in the proton kinetic energies, 
and corrections for variations in the electron and pro- 
ton path lengths through the spectrometers. Pion re- 
jection was performed using a flight-time cut for 7r+s in 
the HRS/j and the Gas Cerenkov for ir~s in the HRS e . 
A sharp, clear, coincidence Time-of-Flight (TOF) peak 
with a FWHM of 1.8 ns resulted (see Fig. 3). High- 
energy correlated protons which punched through the 
HRS/j collimator (<10% of the prompt yield) were re- 
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FIG. 3: Yield spectrum obtained at -Ebeam = 0.843 GeV and 
9pq = +8°, corresponding to p m i ss = 148 MeV/c. Pion rejec- 
tion has been performed, and all timing corrections have been 
applied. The top panel shows a scatterplot of p m i ss versus 
-Bmiss. The dark vertical bands project into the peaks located 
at 12.1 and 18.3 MeV in the bottom panel. These peaks corre- 
spond to protons knocked-out of the lpi/2- and lp3/2-states of 
16 0, respectively. The S m iss resolution was roughly 0.9 MeV 
FWHM, which did not allow for separation of the 2si/2 1^5/2- 
doublet located at -E m i ss = 17.4 MeV from the lp 3 /2-state at 
18.3 MeV. The bump located at roughly 23 MeV is a negative- 
parity doublet which was not investigated. The insert shows 
the corresponding optimized coincidence TOF peak which has 
a FWHM of 1.8 ns. The signal-to-noise ratio was about 8:1 
in these kinematics. 



jected by requiring both spectrometers to independently 
reconstruct the coincidence-event vertex in the vicinity of 
the same water film. The resulting prompt-peak yields 
for each water film were corrected for uncorrelated (ran- 
dom) events present in the peak-time region on a bin- 
by-bin basis as per the method suggested by Owens [62] . 
These per-film yields were then normalized individually. 



B. Normalization 

The relative focal-plane efficiencies for each of the two 
spectrometers were measured independently for each of 
the three water films at every spectrometer excitation 
used in the experiment. By measuring the same single- 
arm cross section at different locations on the spectrome- 
ter focal planes, variations in the relative efficiencies were 
identified. The position variation across the focal plane 



was investigated by systematically shifting the central 
excitation of the spectrometer about the mean momen- 
tum setting in a series of discrete steps such that the full 
momentum acceptance was 'mapped'. A smooth, slowly 
varying dip-region cross section was used instead of a 
single discrete peak for continuous coverage of the focal 
plane. The relative-efficiency profiles were unfolded from 
these data using the program releff [63] by Baghaei. 
For each water film, solid-angle cuts were then applied to 
select the flat regions of the angular acceptance. These 
cuts reduced the spectrometer apertures by roughly 20% 
to about 4.8 msr. Finally, relative-momentum cuts were 
applied to select the flat regions of momentum accep- 
tance. These cuts reduced the spectrometer momentum 
acceptance by roughly 22% to -3.7% < S < 3.3%. The 
resulting acceptance profile of each spectrometer was uni- 
form to within 1%. 

The absolute efficiency at which the two spectrometers 
operated in coincidence mode was given by 



£ ' £p ^coim (1) 

where e e was the single-arm HRS e efficiency, e p was the 
single- arm HRS^ efficiency, and e C om was the coincidence- 
trigger efficiency. The quantity (e p ■ e co i n ) was measured 
at 9 pq = 0° at £ b cam = 0.843 GeV using the ^(e, e) 
reaction. A 0.7 msr collimator was placed in front of 
the HRS e . fn these kinematics, the cone of recoil pro- 
tons fit entirely into the central flat-acceptance region 
of the HRS/j. The number of 1 H(e, e) events where the 
proton was also detected was compared to the number 
of 1 H(e, e) events where the proton was not detected to 
yield a product of efficiencies (e p • e CO i n ) of 98.9%. The 
1.1% effect was due to proton absorption in the water- 
fall target exit windows, spectrometer windows, and the 
first layer of trigger scintillators. Since the central field 
of the HRSft, was held constant throughout the entire ex- 
periment, this measurement was applicable to each of the 
hadron kinematics employed. A similar method was used 
to determine the quantity (e e • e C oin) at each of the three 
HRS e field settings. Instead of a collimator, software 
cuts applied to the recoil protons were used to ensure 
that the cone of scattered electrons fit entirely into the 
central flat-acceptance region of the HRS e . This product 
of efficiencies was >99%. Thus, the coincidence efficiency 
fcoin was firmly established at nearly 100%. A nominal 
systematic uncertainty of ±1.5% was attributed to e. 

The quantity (L ■ e e ), where L is the luminosity (the 
product of the effective target thickness and the number 
of incident electrons) was determined to ±4% by com- 
paring the measured 1 H(e, e) cross section for each film 
at each of the electron kinematics to a parametrization 
established at a similar Q 2 by Simon et al. [64] and Price 
et al. [65] (see Fig. 4). The results reported in this paper 
have all been normalized in this fashion. As a consistency 
check, a direct absolute calculation of (L • e e ) using infor- 
mation from the BCMs, the calibrated thicknesses of the 
water films, and the single-arm HRS e efficiency agrees 
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fit is 1.00 ±0.01 



+ t 



0.843 GeV 1 .643 GeV 



2.442 GeV 



Corrections based on the TS output-to-input ratio 
were applied to the data to account for the acquisition 
deadtime to coincidence events. On average, these cor- 
rections were roughly 20%. An acquisition Monte Carlo 
by Liang [66] was used to cross-check these corrections 
and establish the absolute uncertainty in them at 2%. 

Corrections to the per-film cross-section data for elec- 
tron radiation before and after scattering were calculated 
on a bin-by-bin basis in two ways: first using a version 
of the code RADCOR by Quint [67] modified by Florizone 
[68], and independently, the prescriptions of Borie and 
Dreschel [69] modified by Templon et al. [70] for use 
within the simulation package mceep written by Ulmcr 
[71] . The two approaches agreed to within the statistical 
uncertainty of the data and amounted to <55% of the 
measured cross section for the bound states, and <15% 
of the measured cross section for the continuum. Correc- 
tions for proton radiation at these energies are much less 
than 1% and were not performed. 



time [arb. units] 



FIG. 4: Measured ^(e, e) cross-section data normalized to 
the absolute predictions of a parametrization at similiar Q 2 . 
Statistical error bars are shown. The data shown were taken 
over the course of a three-month run period. The different 
data points for each Sbeam represent different HRSh angular 
settings. 

within uncertainty. 

At every kinematics, a Monte Carlo of the phase-space 
volume subtended by each experimental bin was per- 
formed. For each water foil, No software (e, e'p) events 
were generated, uniformly distributed over the scattered- 
electron and knocked-out proton momenta (p e , p p ) and 
in-plane and out-of-plane angles (</> e , ei cf> p , 6 p ). For 
each of these events, all of the kinematic quantities were 
calculated. The flat-acceptance cuts determined in the 
analysis of the relative focal-plane efficiency data were 
then applied, as were all other cuts that had been per- 
formed on the actual data. The pristine detection volume 
AVb(E miss , p miss , u, Q 2 ) subtended by a bin b(AE miss , 
Ap m i ss , Auj, AQ 2 ) containing Nj, pseudoevents was thus 

AV b (E miss ,p miss ,u, Q 2 ) = 

^ [(Ap e • Aft e ) • (A Pp • Afi p )] , (2) 

where the quantity (Ap e ■ Afi e ) • (Ap p ■ Af2 p ) was the to- 
tal volume sampled over in the Monte Carlo (purposely 
set larger than the experimental acceptance in all dimen- 
sions [139]). The pseudodata were binned exactly as the 
real data, and uniformly on both sides of q. At each 
kinematics, the bin with the largest volume AVJ nax was 
located. Only bins subtending volumes larger than 50% 
of AVJ nax were analyzed further. 



C. Cross section 

The radiativcly corrected average cross section in the 
bin b(AE m i ss , Ap m j ss , Aw, AQ 2 ) was calculated accord- 
ing to 

/ £ \ . 

\ duj dUe dismiss dfl p I b 

Rw O(e,e'p) ( Yb \ ^ 

^coin 

where Yf, was the total number of real events which were 
detected in b(AE miss , Ap miss , Am, AQ 2 ), AV& was the 
phase-space volume, and Rieo( e ,e'p) was a correction ap- 
plied to account for events which radiated in or out of 
AVfc. The average cross section was calculated as a func- 
tion of -E m i ss for a given kinematic setting [140]. Bound- 
state cross-section data for the lp-shell were extracted by 
integrating over the appropriate range in E m i SSl weighting 
with the appropriate Jacobian [141]. Five-fold differen- 
tial cross-section data for QE proton knockout from the 
lp-shell of 16 O are presented in Tables XI and XII. Six- 
fold differential cross-section data for QE proton knock- 
out from 16 O at higher E^ss are presented in Tables XIII 
- XXII. 



D. Asymmetries and response functions 

In the One-Photon Exchange Approximation, the un- 
polarizcd six-fold differential cross section may be ex- 
pressed in terms of four independent response functions 
as (see Refs. [2, 9, 72]) 
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dui d£l e dE m iss d^lp 



where K is a phase-space factor, ciMott is the Mott cross 
section, and the v\ are dimensionless kinematic factors 
[142]. Ideal response functions are not directly measure- 
able because electron distortion does not permit the az- 
imuthal dependences to be separated exactly. The ef- 
fective response functions which are extracted by apply- 
ing Eq. (4) to the data are denoted Rl (longitudinal), 
Rt (transverse) , Rlt (longitudinal-transverse) , and Rtt 
(transverse-transverse). They contain all the information 
which may be extracted from the hadronic system using 
(e,e'p). Note that the Vi depend only on (cj, Q 2 , 8 e ), 
while the response functions depend on (uj, Q 2 , E m i 8S , 

Pmiss) ■ 

The individual contributions of the effective response 
functions may be separated by performing a series of 
cross-section measurements varying Vi and/or </>, but 
keeping q and w constant [143]. In the case where the 
proton is knocked-out of the nucleus in a direction par- 
allel to q ('parallel' kinematics), the interference terms 
Rlt an d Rtt vanish, and a Rosenbluth separation [73] 
may be performed to separate Rl and Rt- In the case 
where the proton is knocked-out of the nucleus in the 
scattering plane with a finite angle 9 pq with respect to q 
('quasiperpendicular' kinematics), the asymmetry Alt 
and the interference Rlt may be separated by perform- 
ing symmetric cross-section measurements on either side 
oiq{<t> = 0° and <j> = 180°). The contribution of Rtt can- 
not be separated from that of Rl with only in-plane mea- 
surements; however, by combining the two techniques, an 
interesting combination of response functions Rt, Rlt, 
and Rl+tt [144] may be extracted. 

For these data, effective response-function separations 
were performed where the phase-space overlap between 
kinematics permitted. For these separations, bins were 
selected only if their phase-space volumes AVf, were all 
simultaneously 50% of AV^ax- Separated effective re- 
sponse functions for QE proton knockout from the lp- 
shell of 16 are presented in Tables XXIII, XXIV, and 
XXV. Separated effective response functions for QE pro- 
ton knockout from the 16 continuum are presented in 
Tables XXVI, XXVII, and XXVIII. 



E. Systematic uncertainties 

The systematic uncertainties in the cross-section mea- 
surements were classified into two categories - kinematic- 
dependent uncertainties and scale uncertainties. For a 
complete discussion of how these uncertainties were eval- 
uated, the interested reader is directed to a report by 
Fissum and Ulmer [74]. For the sake of completeness, 



vtRt + vltRlt cos(0) + vttRtt cos(20)] , (4) 
I 



TABLE III: Kinematic-dependent systematic uncertainties 
folded into the mceep simulation series. 



Quantity 


description 


5 


E b cam 


beam energy 


1.6 x 1(T 3 


^bcam 


in-plane beam angle 


ignored" 


^bcam 


out-of-plane beam angle 


2.0 mrad 


Pe 


scattered electron momentum 


1.5 x 1(T 3 


</>,: 


in-plane scattered electron angle 


0.3 mrad 


0,: 


out-of-plane scattered electron angle 


2.0 mrad 


Pp 


proton momentum 


1.5 x io~ 3 


4> P 


in-plane proton angle 


0.3 mrad 


Op 


out-of-plane proton angle 


2.0 mrad 



"As previously mentioned, the angle of incidence of the electron 
beam was determined using a pair of BPMs located upstream of 
the target (see Fig. 2). The BPM readback was calibrated by 
comparing the location of survey fiducials along the beamline to 
the Hall A survey fiducials. Thus, in principle, uncertainty in the 
knowledge of the incident electron-beam angle should be included 
in this analysis. However, the simultaneous measurement of the 
kinematically overdctermined 1 H(e, ep) reaction allowed for a cali- 
bration of the absolute kinematics, and thus an elimination of this 
uncertainty. That is, the direction of the beam defined the axis 
relative to which all angles were measured via 1 H(e, ep). 



a subset of the aforementioned information is presented 
here. 

In a series of simulations performed after the exper- 
iment, MCEEP was used to investigate the intrinsic be- 
havior of the cross-section data when constituent kine- 
matic parameters were varied over the appropriate ex- 
perimentally determined ranges presented in Table III. 
Based on the experimental data, the high-i? m i SS region 
was modelled as the superposition of a peak-like lS]/2- 
state on a flat continuum. Contributions to the system- 
atic uncertainty from this flat continuum were taken to 
be small, leaving only those from the ls^ -8 ! 8 -! - The 
16 0(e, e'p) simulations incorporated as physics input the 
bound-nucleon rdwia calculations detailed in Section 
IV A, which were based on the experimental lp-shcll 
data. 

For each kinematics, the central water foil was con- 
sidered, and 1M events were generated. In evaluating 
the simulation results, the exact cuts applied in the ac- 
tual data analyses were applied to the pseudo-data, and 
the cross section was evaluated for the identical p m i 8s bins 
used to present the results. The experimental constraints 
to the kinematic-dependent observables afforded by the 
overdetermined 1 H(e,ep) reaction were exploited to cali- 
brate and constrain the experimental setup. The in-plane 
electron and proton angles 4> e and 4> p were chosen as in- 
dependent parameters. When a known shift in (f> e was 
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made, <p p was held constant and the complementary vari- 
ables -Ebcam: Pe, and p p were varied as required by the 
constraints enforced by the 1 H(e, ep) reaction. Similarly, 
when a known shift in (f> p was made, <f> e was held con- 
stant and the complementary variables -Ebcam, Pe, and 
p p were varied as appropriate. The overall constrained 
uncertainty was taken to be the quadratic sum of the two 
contributions. 

The global convergence of the uncertainty estimate was 
examined for certain extreme kinematics, where 10M- 
event simulations (which demonstrated the same behav- 
ior) were performed. The behavior of the uncertainty 
as a function of p m i ss was also investigated by examin- 
ing the uncertainty in the momentum bins adjacent to 
the reported momentum bin in exactly the same fashion. 
The kinematically induced systematic uncertainty in the 
16 0(e,e'p) cross-section data was determined to be de- 
pendent upon pmiss, with an average value of 1.4%. The 
corresponding uncertainties in the ^(eje) cross-section 
data were determined to be negligible. 

The scale systematic uncertainties which affect each 
of the cross-section measurements are presented in Ta- 
ble IV. As previously mentioned, the 16 0(e,e'p) cross- 
section results reported in this paper have been nor- 
malized by comparing simultaneously measured 1 H(e, e) 
cross-section data to a parametrization established at a 
similar Q 2 . Thus, the first seven listed uncertainties sim- 
ply divide out of the quotient, such that only the subse- 
quent uncertainties affect the results. The average sys- 
tematic uncertainty associated with a lp-shell cross sec- 
tion was 5.6%, while that for the continuum was 5.9%. 
The small difference was due to contamination of the 
high-i^miss data by collimator punch-through events. 

The quality of these data in terms of their associated 
systematic uncertainties was clearly demonstrated by the 
results obtained for the effective response-function sepa- 
rations. In Fig. 5, cross-section data for the lp-shell 
measured in parallel kinematics at three different beam 
energies are shown as a function of the separation lever 
arm vt/vl- The values of the effective response func- 
tions Rl (offset) and Rt (slope) were extracted from the 
fitted line. The extremely linear trend in the data indi- 
cated that the magnitude of the systematic uncertainties 
was small, and that statistical uncertainties dominated. 
This is not simply a test of the One-Photon Exchange 
Approximation employed in the data analysis as it has 
been demonstrated by Traini et al. [75] and Udias [76] 
that the linear behavior of the Roscnbluth plot persists 
even after Coulomb distortion is included. 

Given the applicability of the One-Photon Exchange 
Approximation at these energies, the quality of the data 
was also demonstrated by the results extracted from iden- 
tical measurements which were performed in different 
electron kinematics. The asymmetries Alt and effective 
response functions Rlt for QE proton knockout were ex- 
tracted for both £; boam = 1.643 GcV and 2.442 GcV for 
9 pq = ±8° (p miss = 148 MeV/c). They agree within the 
statistical uncertainty. Table XXV presents the results 




FIG. 5: Cross-section data for the removal of protons from 
the lp-shell of 16 O measured in parallel kinematics at three 
different beam energies as a function of the separation lever 
arm vt/vl- The data points correspond to beam energies of 
2.442, 1.643, and 0.843 GeV from left to right. The effec- 
tive response functions Rl (offset) and Rt (slope) have been 
extracted from the fitted line. The uncertainties shown are 
statistical only. The extremely linear behavior of the data 
(which persists even after corrections for Coulomb distortion 
are applied) indicates that the statistical uncertainties were 
dominant (see Section III E for a complete discussion) . 

at both beam energies for lp-shell knockout for < Q 2 > 
= 0.800 (GcV/c) 2 , < u> > = 436 MeV, and < T p > = 
427 MeV, while Fig. 6 shows the results for 25 < E miss 
< 60 MeV. 



IV. THEORETICAL OVERVIEW 

In the following subsections, overviews of Relativistic 
Distorted- Wave Impulse Approximation (rdwia), Rela- 
tivistic Optical-Model Eikonal Approximation (romea), 
and Relativistic Multiple-Scattering Glauber Approxi- 
mation (rmsga) calculations are presented. 

A. RDWIA 

Reviews of work on proton electromagnetic knockout 
using essentially nonrelativistic approaches may be found 
in Refs. [1-3]. As previously mentioned, the Relativis- 
tic Distorted- Wave Impulse Approximation (rdwia) was 
pioneered by Picklesimer, Van Orden, and Wallace [7- 
9] and subsequently developed in more detail by several 
groups (see Refs. [16, 21, 77-81]). In Section IV A 1, 
the RDWIA formalism for direct knockout based upon 
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TABLE IV: Summary of the scale systematic uncertainties contributing to the cross-section data. The first seven entries do 
not contribute to the systematic uncertainties in the reported cross-section data as they contribute equally to the 1 H(e,e) 
cross-section data to which the 16 0(e,e'p) are normalized. 



Quantity 


description 


8(%) 


??DAQ 


data acquisition deadtime correction 


2.0 




electronics deadtime correction 


<1.0 


pt' 


effective target thickness 


2.5 


N c 


number of incident electrons 


2.0 


e e 


electron detection efficiency 


1.0 




HRS e solid angle 


2.0 


£e ' ' Ccoin 


product of electron, proton, and coincidence efficiencies 


1.5 


L ■ e e 


obtained from a form-factor parametrization of 1 H(e,e) 


4.0 


■Kl60(e,e'j>) 


radiative correction to the le O(e,e'p) data 


2.0 


-R 1 H(e,e) 


radiative correction to the 1 H(e, e) data 


2.0 


Cp ' Ccoin 


product of proton and coincidence efficiencies 


<1.0 


AQ p a 


HRSh solid angle 


2.0 


punchthrough c 


protons which punched through the HRSh collimator 


2.0 



"The systematic uncertainties in the solid angles AQ e and AQ P were quantified by studying sieve-slit collimator optics data at each of 
the spectrometer central momenta employed. The angular locations of each of the reconstructed peaks corresponding to the 7x7 lattice 
of holes in the sieve-slit plate were compared to the locations predicted by spectrometer surveys, and the overall uncertainy was taken to 
be the quadratic sum of the individual uncertainties. 

'At first glance, it may be surprising to note that the uncertainty due to the radiative correction to the data is included as a scale 
uncertainty. In general, the radiative correction is strongly dependent on kinematics. However, the lp-shell data analysis, and for that 
matter any bound-state data analysis, involves -E m iss cuts. These cuts to a large extent remove the strong kinematic dependence of the 
radiative correction, since only relatively small photon energies arc involved. In order to compensate for any remaining weak kinematic 
dependence, the uncertainty due to the radiative correction was slightly overestimated. 

c High -E m ; ss data only. 
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FIG. 6: Rlt for 8 pq = ±8° (p miss = 145 MeV/c) as a function 
of E miBS for £ bcam = 1.643 GeV and 2.442 GeV. Statistical 
uncertainties only are shown. The statistical agreement over 
a broad range of -Emiss emphasizes the systematic precision of 
the measurement (see Section HIE for a complete discussion). 
Note that the averages of these Rlt values are presented as 
the pmiss = 145 MeV/c data in Fig. 22 and Table XXVIII. 



a singlc-nucleon operator is outlined in sufficient detail 
that the most important differences with respect to non- 
relativistic dwia may be identified. In Section IV A 2, 
a direct numerical comparison between two different im- 
plementations of rdwia is presented. 



1. Formalism 

The five-fold differential cross section for the exclusive 
A(e, e'N)B reaction leading to a discrete final state takes 
the form (see Ref. [2]) 



d 5 o- 



= K 



N 



et Q 



where 



K = K 



PnE_ 



N 



(2nf 



(5) 



(6) 



is a phase-space factor, fc, = (ej, fc,) and kt — (ef,kf) 
are the initial and final electron momenta, pa = (Ea,Pa) 
and pb = (Eb,Pb) are the initial and final target mo- 
menta, pn = {En,Pn) is the ejected-nucleon momen- 
tum, q = ki — kf = (u),q) is the momentum transfer 
carried by the virtual photon, Q 2 = —q^.q^ = q 2 — 
the photon virtuality, and 



U) is 



1 - 



vn ■ vb 
vn ■ v N 



(7) 
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(with i'n = pn/En) is a recoil factor which adjusts the 
nuclear phase space for the missing-energy constraint. In 
the One-Photon Exchange Approximation, the invariant 
electroexcitation matrix element is represented by the 
contraction of electron and nuclear response tensors of 
the form 



fit/ 



(8) 
(9) 



where j M is the electron current, is a matrix element 
of the nuclear electromagnetic current, and the angled 
brackets denote averages over initial states and sums over 
final states. 

The reduced cross section is given by 



d 5 a 



where 



Crcd 



G eN 



-/Ka 



eN , 



e/ a 



PWIA 



(10) 



(11) 



is the elementary cross section for electron scattering 
from a moving free nucleoli in the Plane- Wave Impulse 
Approximation (pwia). The PWIA response tensor is 
computed for a free nucleon in the final state, and is 
given by 



Kwia = ^TraccJ^t , 



where 



u(p f ,Sf)T^u(pi,Si 



(12) 



(13) 



is the single- nucleon current between free spinors normal- 
ized to unit flux. The initial momentum (pi = pf — q e s) 
is obtained from the final ejectile momentum (pf) and 
the effective momentum transfer (g c ff) in the laboratory 
frame, and the initial energy is placed on shell. The 
effective momentum transfer accounts for electron accel- 
eration in the nuclear Coulomb field and is discussed fur- 
ther later in this Section. 

In the nonrclativistic PWIA limit, er roc i reduces to the 
bound-nuclcon momentum distribution, and the cross 
section given in Eq. (5) may be expressed as the product 
of the phase-space factor K, the elementary cross section 
a e N, and the momentum distribution. This is usually 
referred to 'factorization'. Factorization is not strictly 
valid relativistically because the binding potential alters 
the relationship between lower and upper components of 
a Dirac wave function - see Ref . [82] . 

In this Section, it is assumed that the nuclear current 
is represented by a one-body operator, such that 



J" 



d 3 r exp(^-r)(* ( ~ ) (p,r)|r Al |< ? !)(r)) , (14) 



where <fi is the nuclear overlap for singlc-nucleon knockout 
(often described as the bound-nucleon wave function), 
is the Dirac adjoint of the time-reversed distorted 
wave, p is the relative momentum, and 



Eb 
W ' 



(15) 



is the recoil-corrected momentum transfer in the barycen- 
tric frame. Here (u>, q ) and Eb are the momentum trans- 
fer and the total energy of the residual nucleus in the lab- 



oratory frame respectively, and W = y/ {ttia + u) 2 q 2 
is the invariant mass. 

De Forest [83] and Chinn and Picklesimer [84] have 
demonstrated that the electromagnetic vertex function 
for a free nucleon can be represented by any of three 
Gordon-equivalent operators 



pp. 

^ 1 (Pf,p l )^Y l G M (Q 2 ) - ^F 2 (Q 2 ) 

r%(p f ,p l )=Y l F 1 (Q 2 )+z<j^^-F 2 (Q 



(16a) 
(16b) 



PM n 

r$(p f ,Pi) = ^ F i(Q 2 ) + ^7^g m (Q 2 ) (i6c) 

where P = (Ef + E i7 pf + pi). Note the correspon- 
dence with Eq. (17) below. Although T 2 is arguably 
the most fundamental because it is defined in terms of 
the Dirac and Pauli form factors F\ and F 2 , Ti is often 
used because the matrix elements are easier to evaluate. 
T3 is rarely used but no less fundamental. In all calcula- 
tions presented here, the momenta in the vertex functions 
are evaluated using asymptotic laboratory kinematics in- 
stead of differential operators. 

Unfortunately, as bound nucleons are not on shell, an 
off-shell extrapolation (for which no rigorous justifica- 
tion exists) is required. The de Forest prescription is 
employed, in which the energies of both the initial and 
the final nucleons are placed on shell based upon effec- 
tive momenta, and the energy transfer is replaced by the 
difference between on-shell nucleon energies in the oper- 
ator. Note that the form factors arc still evaluated at the 
Q 2 determined from the electron-scattering kinematics. 
In this manner, three prescriptions 



= rG M {Q 2 )~—F 2 {Q 2 ) 
2m 

f£ = ^ Fl {Q 2 )+ia^^-F 2 (Q 2 ) 

n = —F l (Q 2 ) +l a^^G M {Q 2 ), (17c) 
2m 2m 



(17a) 
(17b) 



arc obtained, where 



q = (E f - Ei,q) , 
P = (E f + Ei,2p f - q) , 



and where Ei — \Jm 2 N + (pf — q) 2 is placed on shell 
based upon the externally observable momenta p f and q 
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evaluated in the laboratory frame. When electron distor- 
tion is included, the local momentum transfer q — > q ff 
is interpreted as the effective momentum transfer with 
Coulomb distortion. These operators are commonly 
named CCl, CC2, and CC3, and are no longer equiva- 
lent when the nucleons are off-shell. Furthermore, the ef- 
fects of possible density dependence in the nucleon form 
factors can be evaluated by applying the Local Density 
Approximation (LDA) to Eq. (17) - see Refs. [80, 84]. 

The overlap function is represented Dirac spinor 
of the form 



(f> Km {r) 



f K (r)y Km (f) 



(18) 



where 



y Km {f) = ^2 / 1 v - m s \j m \ Yi u (f)x,n s (19) 



is the spin spherical harmonic and where the orbital and 
total angular momenta are respectively given by 



i = s k k 



l 



(20a) 
(20b) 



with S K = sign(«;). The functions f K and g K satisfy the 
usual coupled linear differential equations - see for ex- 
ample Rcf. [85]. The corresponding momentum wave 
function 

4> Km (Pm ) = J d 3 r exp (-ip m ■ r )<f) Km (r ) (21) 
then takes the form 

fK(Pm)yKm(Pm) 



6 (n )=Am~ t l Jn\l>m)J'Km\t>m) /oo\ 

<p Km {Pm) 4m [ .s K g- K {p m )y- Km (p m ) ' ' [22) 



where 



f K (Pm) = dr r 2 j e (p m r)f K (r) (23a) 



and where in the PWIA, the initial momentum p m would 
equal the experimental missing momentum p m i ss . Thus, 
the momentum distribution 



P{Pm) 



1 



2tt 2 

is obtained, normalized to 



\f K {pm)\ 2 + \g K (p m )\ 



4?r J dpp 2 m p{p m ) = 1 



(24) 



(25) 



for unit occupancy. 



Similarly, let 



#(+)(p,r) 



2E 



CM 



(26) 



represent a wave function of the N+B system with an in- 
coming Coulomb wave and outgoing spherical waves open 
in all channels. Specific details regarding the boundary 
conditions may be found in Refs. [86-88]. 

The Madrid rdwia calculations [16] employ a partial- 
wave expansion of the first-order Dirac equation, leading 
to a pair of coupled first-order differential equations. Al- 
ternatively, the LEA code [89] by Kelly uses the Numerov 
algorithm to solve a single second-order differential equa- 
tion that emerges from an equivalent Schrodinger equa- 
tion of the form 



k 2 



2(i (U c + U LS L ■ tr)] £ = 



(27) 



where k is the relativistic wave number, /i is the reduced 
energy, and 



E 
1 



s 2 -v 2 

2E 



U 



D 



u 



LS 



D = 1 



2r 2 D dr 

_iy_ 

s-v 



1 



E 



(28a) 

(28b) 

(28c) 
(28d) 



S and V [145] are respectively the scalar and vector po- 
tential terms of the original four-component Dirac equa- 
tion (sec Rcf [2]). D(r) is known as the Darwin nonlocal- 
ity factor and U c and U LS are the central and spin-orbit 
potentials. The Darwin potential U D is generally quite 
small. The upper and lower components of the Dirac 
wave function are then obtained using 





c 



D l/2£ 

er • p ip 



E 



S-V 



(29a) 
(29b) 



This method is known as direct Pauli reduction [18, 79]. 
A very similar approach is also employed by Meucci et al. 
[81]. A somewhat similar approach based on the Eikonal 
Approximation (see the discussion of the ROMEA calcu- 
lations in Section IV B) has been employed by Radici et 
al. [90, 91]. 

For our purposes, the two most important differences 
between relativistic and nonrclativistic dwia calculations 
are the suppression of the interior wave function by the 
Darwin factor in Eq. (29a), and the dynamical enhance- 
ment of the lower components of the Dirac spinor (also 
known as 'spinor distortion') by the strong Dirac scalar 
and vector potentials in Eq. (29b). 

As demonstrated in Refs. [18, 92, 93], the Darwin 
factor tends to increase the normalization factors de- 
duced using an rdwia analysis. Distortion of the bound- 
nucleon spinor destroys factorization and at large p m i ss 
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produces important oscillatory signatures in the interfer- 
ence response functions, Alt, an d recoil polarization - 
see Rcfs. [19 21, 88, 94]. The effect of spinor distortion 
within the Effective Momentum Approximation (EMA) 
has been studied by Kelly [88]. The LEA code has sub- 
sequently been upgraded to evaluate Eq. (29) without 
applying the EMA. These two methods for constructing 
the ejectile distorted waves should be equivalent. The 
predictions of the LEA and the Madrid codes given iden- 
tical input are compared in Section IV A 2. 

The approximations made by dwia violate current 
conservation and introduce gauge ambiguities. The most 
common prescriptions 



quasiperpendicular 



parallel 



J, 
Jo 



-Jo 

q 



J, 
q 



J_q 

Q 2 



J, 



(30a) 
(30b) 
(30c) 



correspond to Coulomb, Landau, and Weyl gauges, re- 
spectively. Typically, Gordon ambiguities and sensitivity 
to details of the off-shell extrapolation are largest in the 
Weyl gauge. Although there is no fundamental prefer- 
ence for any of these prescriptions, it appears that the 
data are in general least supportive of the Weyl gauge. 
Further, the CCl operator is the most sensitive to spinor 
distortion while the CC3 operator is the least. The inter- 
mediate CC2 is chosen most often for rdwia. 

Besides the interaction in the final state of the outgoing 
proton, in any realistic calculation with finite nuclei, the 
effect of the distortion of the electron wave function must 
be taken into account. For relatively light nuclei and 
large kinetic energies, the EMA for electron distortion 
(the (7 G ff Approximation) is sufficient - see Rcfs. [95, 96]. 
In this approach, the electron current is approximated by 



kikj 
hi k f 



(31) 



where q c g = ki — kf is the effective momentum transfer 
based upon the effective wave numbers 



k = k + f z ^k 

Hz 



(32) 



with fz ps 1.5 and Rz ps 1.2A 1/3 . For all the RDWIA cal- 
culations subsequently presented in this paper that are 
compared directly with data, this q a R Approximation has 
been used to account for electron Coulomb distortion. 
Only the RPWIA and RDWiA(uo P t =0 ) comparison calcula- 
tions shown in Fig. 7, the baseline rdwia comparison 
calculations shown in Figs. 8 and 9, and the rdwia and 
RMSGA comparison calculations shown in Fig. 10 omit 
the effect of electron Coulomb distortion. This is equiv- 
alent to setting fz = in Eq. (32). 
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FIG. 7: Comparisons between rpwia and RDWlA( U o P t =0 j cal- 
culations for the removal of protons from the lpx/2-state 
of 16 O as a function of p m i ss for -Ebcam = 2.442 GeV. In 
the upper panels, the solid curves represent the reduced 
cross sections for both the RDWlA( U o P t =0 - ) and the rpwia 
calculations (see text for details). The dashed curves cor- 
respond to the momentum distributions. In the lower pan- 
els, RDWlA(uopt =0 ) / RPWIA reduced cross-section ratios are 
shown. Agreement to much better than 1% is obtained for 
both kinematics over the entire p m iss range. 



2. Tests 

Fig. 7 illustrates a comparison between RPWIA and 
RDWlA(uopt =0 ) calculations made by Kelly using LEA for 
the removal of protons from the lpi/2-state of 16 as a 
function of p m i ss for both quasiperpendicular and parallel 

(U o p'=0) 



kinematics for E^ c 



2.442 GeV. The rdwia 



calculations employed a partial-wave expansion of the 
second-order Dirac equation with optical potentials nul- 
lified and the target mass artificially set to 16001m to 
minimize recoil corrections and frame ambiguities. The 
rpwia calculations (see Ref. [5]) are based upon the 
Fourier transforms of the upper and lower components of 
the overlap function; that is, no partial-wave expansion 
is involved. In the upper panels, the solid curves repre- 
sent the reduced cross section for both the RDWiA(u°p'=o) 
and RPWIA calculations as the differences are indistin- 
guishable on this scale. The dashed curves show the mo- 
mentum distributions. In the lower panels, the ratios be- 
tween RDWiA(Tjo P t =0 ) and RPWIA reduced cross sections 
are shown. With suitable choices for step size and maxi- 
mum £ (here 0.05 fm and 80), agreement to much better 
than 1% over the entire range of missing momentum is 
obtained, verifying the accuracy of LEA for plane waves. 
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TABLE V: A summary of the basic rdwia options which 
served as input to the 'baseline' comparison calculations of 
the Madrid Group and Kelly (LEA). Results are shown in Figs. 
8 and 9. 



Input Parameter 



Option 



bound-nucleon wave function 
Optical Model 
nucleon spinor distortion 
electron distortion 
current operator 
nucleon form factors 
gauge 



nlsh - Sharma et al. [98] 
edai-O - Cooper et al. [99] 
relativistic 
none 
CC2 
dipole 
Coulomb 



Similar results arc obtained with the Madrid code. 

The similarity between the reduced cross sections and 
the momentum distributions demonstrates that the vio- 
lation of factorization produced by the distortion of the 
bound-state spinor is mild, but tends to increase with 
Pmiss- Nevertheless, observables such as Alt that are 
sensitive to the interference between the lower and up- 
per components are more strongly affected by violation 
of factorization. 

Figs. 8 and 9 compare calculations for the removal of 
protons from the lp-shell of 16 O as a function of p m i ss 
for E'bcam = 2.442 GeV. These predictions were made 
by Kelly [97] using the LEA and the Madrid codes with 
identical input options, and are hereafter described as 
'baseline' calculations. The baseline options are summa- 
rized in Table V and were chosen to provide the most 
rigorous test of the codes rather than to be the optimal 
physics choices. 

Fig. 8 demonstrates that baseline cross-section calcula- 
tions agree to better than 2% for p m i ss < 250 MeV/c, but 
that the differences increase to about 10% by about 400 
MeV/c. Nevertheless, Fig. 9 shows that excellent agree- 
ment is obtained for Alt over this entire range of p m i ss , 
with only a very small observable shift. The agreement 
of the strong oscillations in Alt for Pmiss ~ 300 MeV/c 
predicted by both methods demonstrates that they are 
equivalent with respect to spinor distortion. The small 
differences in the cross section for large p m i ss appear to be 
independent of the input choices and probably arise from 
numerical errors in the integration of differential equa- 
tions (perhaps due to initial conditions), but the origin 
has not yet been identified. Regardless, it is remarkable 
to achieve this level of agreement between two indepen- 
dent codes under conditions in which the cross section 
spans three orders of magnitude. 



B. ROMEA / RMSGA 

In this subsection, an alternate relativistic model de- 
veloped by the Ghent Group [100-103] for A(e, e'N)B 
processes is presented. With respect to the construction 
of the bound-nucleon wave functions and the nuclear- 
current operator, an approach similar to standard rdwia 
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FIG. 8: Comparison baseline rdwia calculations by the 
Madrid Group and Kelly (LEA) for the removal of protons 
from the lp-shell of 16 as a function of p m i S s for -Ebcam = 
2.442 GeV. For the purposes of this comparison, the input 
into both calculations was identical (see Table V). Overall 
agreement is very good, and agreement is excellent for —250 
< Pmiss < 250 MeV/c. 
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FIG. 9: Comparison baseline rdwia calculations for the Alt 
asymmetry by the Madrid Group and Kelly (LEA) for the re- 
moval of protons from the lp-shell of 16 O as a function of p m i ss 
for -Ebeam = 2.442 GeV. For the purposes of this comparison, 
the input into both calculations was identical (see Table V). 
Overall agreement is excellent over the entire p m i ss range. 
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is followed. The major differences lie in the construction 
of the scattering wave function. The approach presented 
here adopts the relativistic Eikonal Approximation (EA) 
to determine the scattering wave functions and may be 
used in conjunction with either the Optical Model or the 
multiple-scattering Glauber frameworks for dealing with 
the FSI. 



1. Formalism 

The EA belongs to the class of semi-classical approxi- 
mations which arc meant to become 'exact' in the limit 
of small de Broglic (db) wavelengths. Xdb <C a, where a 
is the typical range of the potential in which the particle 
is moving. For a particle moving in a relativistic (opti- 
cal) potential consisting of scalar and vector terms, the 
scattering wave function takes on the EA form 



^ A {r,r2,r 3 , ...r A ) 



O 



1 



x^ B (r 2l ...r A ) , (35) 



where is the wave function characterizing the state 
in which the B nucleus is created. In the above expres- 
sion, the subsequent elastic or 'mildly inelastic' collisions 
which the ejectile undergoes with 'frozen' spectator nu- 
cleons are implemented through the introduction of the 
operator 



d(r,r 2 ,r 3 ,...r A ) = l[[l-T{p,b-b j )e(z 
where the profile function for pN scattering is 



°0)\ J 



ip F (r) 



E+m+S-V 



IT ■ p 



ipr iS(r) 



(33) 



This wave function differs from a relativistic plane wave 
in two respects: first, there is a dynamical relativistic ef- 
fect from the scalar (S) and vector (V) potentials which 
enhances the contribution from the lower components; 
and second, the wave function contains an eikonal phase 
which is determined by integrating the central (U ) and 
spin-orbit (U L ) terms of the distorting potentials along 
the (asymptotic) trajectory of the escaping particle. In 
practice, this amounts to numerically calculating the in- 
tegral (r = (b, z)) 



iS(b,z) 



' l K 1 dz 



U c (b,z') 



+U LS {b, z')[cr ■ (b X K) - iKz'} 



(34) 



where K = \{p + q)- 

Within the romea calculation, the eikonal phase given 
by Eq. (34) is computed from the relativistic optical 
potentials as they are derived from global fits to elastic 
proton-nucleus scattering data. It is worth stressing that 
the sole difference between the romea and the rdwia 
models is the use of the EA to compute the scattering 
wave functions. 

For proton lab momenta exceeding 1 GeV/c, the 
highly inelastic nature of the elementary nuclcon-nuclcon 
(NN) scattering process makes the use of a potential 
method for describing FSI effects somewhat artificial. In 
this high-energy regime, an alternate description of FSI 
processes is provided by the Glauber Multiple-Scattering 
Theory. A relativistic and unfactorized formulation of 
this theory has been developed by the Ghent Group 
[102, 103]. In this framework, the ^4-body wave function 
in the final state reads 



r(p, 6) = 



'pN 



(1 - ie pN ) 



4^/3 2 



cxp 



pN 



In practice, for the lab momentum of a given ejectile, 
the following input is required: the total proton-proton 
and proton- neutron cross section o~ P °^, the slope param- 
eters PpN, and the ratio of the real-to- imaginary scat- 
tering amplitude c p n. The parameters c*^, (3 p n, and 
t P N are obtained through interpolation of the data base 
made available by the Particle Data Group [104]. The 
A(e, e'N)B results obtained with a scattering state of the 
form of Eq. (35) are referred to as RMSGA calculations. It 
is worth stressing that in contrast to the rdwia and the 
romea models, all parameters entering the calculation 
of the scattering states in RMSGA are directly obtained 
from the elementary proton-proton and proton-neutron 
scattering data. Thus, the scattering states are not sub- 
ject to the SV effects discussed in Section IV A 1, which 
typically arise when relativistic potentials are employed. 
However, the SV effects are included for the bound-state 
wave function. 

Note that for the kinematics of the 16 0(e,e'p) experi- 
ment presented in this paper, the de Broglie wavelength 
of the ejected proton is Xdb ~ 1-3 fm, and thus both 
the optical potential and the Glauber frameworks may 
be applicable. Indeed, for T p w 0.433 GcV, various sets 
of relativistic optical potentials are readily available and 
Xdb appears sufficiently small for the approximations en- 
tering the Glauber framework to be justifiable - see Ref. 
[102]. 



2. Tests 

As previously mentioned, the (rdwia and romea) and 
RMSGA frameworks are substantially different in the way 
they address FSI. While the rdwia and ROMEA mod- 
els are both essentially one-body approaches in which all 
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TABLE VI: A summary of the basic options which served as 
input to the comparison between the rdwia calculations and 
the 'bare' rmsga (no MEC nor IC) calculations of the Ghent 
Group. Results are shown in Fig. 10. 



Input Parameter 



Option 



bound-nucleon wave function 
Optical Model 
nucleon spinor distortion 
electron distortion 
current operator 
nucleon form factors 
gauge 



Furnstahl et al. [106] 

EDAI-O 

relativistic 
none 
CC2 
dipole 
Coulomb 



FSI effects are implemented through effective potentials, 
the RMSGA framework is a full-fledged, multi-nucleon 
scattering model based on the EA and the concept of 
frozen spectators. As such, when formulated in an un- 
factorized and relativistic framework, Glauber calcula- 
tions are numerically involved and the process of com- 
puting the scattering state and the transition matrix el- 
ements involves numerical methods which are different 
from those adopted in rdwia frameworks. For exam- 
ple, for A(e, e'N)B calculations in the romea and the 
RMSGA, partial-wave expansions are simply not a viable 
option. 

The testing of the mutual consistency of the rdwia 
and 'bare' RMSGA (no MEC nor IC) calculations began 
by considering the special case of vanishing FSI. In this 
limit, where all the Glauber phases are nullified in RMSGA 
and (rdwia — > rpwia), the two calculations were deter- 
mined to reproduce one another to 4% over the entire 
Pmiss range, thereby establishing the validity of the nu- 
merics. The Glauber phases were then enabled. The ba- 
sic options which then served as input to the comparison 
between the rdwia calculations and the RMSGA calcu- 
lations of the Ghent Group [105] are presented in Table 
VI. 

Fig. 10 shows the ratio of the bare RMSGA calculations 
of the Ghent Group together with rdwia calculations 
for the removal of protons from the lp-shell of 16 O as 
a function of p m i ss for -E^am = 2.442 GeV. Apart from 
the treatment of FSI, all other ingredients to the calcu- 
lations are identical (see Table VI). For p m i ss below the 
Fermi momentum, the variation between the predictions 
of the two approaches is at most 25%, with the rdwia 
approach predicting a smaller cross section (stronger ab- 
sorptive effects) than the RMSGA model. Not surprisingly, 
at larger p m ; ss (correspondingly larger polar angles), the 
differences between the two approaches grow. 



V. RESULTS FOR Q 2 « 0.8 (GeV/c) 5 



The data were interpreted in subsets corresponding to 
the lp-shell and to the ls^-statc and continuum. The 
interested reader is directed to the works of Gao et al. 
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FIG. 10: rdwia calculations compared to 'bare' (no MEC 
nor IC) RMSGA calculations by the Ghent Group for the re- 
moval of protons from the lp-shell of 16 O as a function of p m i aa 
for .Ebeam = 2.442 GeV. Both calculations employ the input 
presented in Table VI. Apart from the treatment of FSI, all 
ingredients are identical. 



[107] and Liyanage et al. [108], where these results have 
been briefly highlighted. Note that when data are pre- 
sented in the following discussion, statistical uncertain- 
ties only are shown. A complete archive of the data, 
including systematic uncertainties, is presented in Ap- 
pendix A. 



A. lp-shell knockout 

1. Sensitivity to rdwia variations 

The consistency of the normalization factors S a sug- 
gested by the lp-shell data for p m i ss < 350 MeV/c ob- 
tained in this measurement at 2.442 GeV (see Table XII) 
was examined within the rdwia framework in a detailed 
study by the Madrid Group [109]. The study involved 
systematically varying a wide range of inputs to the rd- 
wia calculations, and then performing least-squares fits 
of the predictions to the cross-section data. The results 
of the study are presented in Table VII. 

Three basic approaches were considered: the fully rel- 
ativistic approach, the projected approach of Udfas et 
al. [20, 21], and the EMA-noSV approach of Kelly 
[4, 88]. All three approaches included the effects of elec- 
tron distortion. While the fully relativistic approach 
involved solving the Dirac equation directly in config- 
uration space, the projected approach included only the 
positive-energy components, and as a result, most (but 
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TABLE VII: Normalization factors derived from the 2.442 GeV lp-shell cross-section data of Table XII using the CCl and CC2 
current operators. The first term in each column is for the lp^-state, while the second term is for the lp 3 / 2 -state. 
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not all) of the spinor distortion was removed from the 
wave functions. Within the EMA-noSV approach, a rel- 
ativized Schrodinger equation was solved using the EMA, 
and all of the spinor distortion was removed. This made 
the calculation similar to a factorized calculation, al- 
though spin-orbit effects in the initial and final states 
(which cause small deviations from the factorized results) 
are included in EMA-noSV. 

The current operator was changed between CCl and 
CC2. Three bound-nuclcon wave functions (see Fig. 11) 
derived from relativistic Lagrangians were considered: HS 
by Horowitz and Serot [110, 111], nlsh by Sharma et 
al. [98], and nlsh-p by Udfas et al. [112] (which re- 
sulted from a Lagrangian fine-tuned to reproduce the 
Leuschner et al. data). Note that both the NLSH and 
NLSH-P wave functions predict binding energies, single- 
particle energies, and a charge radius for 16 which are 
all in good agreement with the data. The gauge prescrip- 
tion was changed between Coulomb, Weyl, and Landau. 
The nucleon distortion was evaluated using three purely 
phenomenological SV optical potentials (edai-O, EDADl, 
and edad2) by Cooper et al. [99], as well as mrw by Mc- 
Neil et al. [113] and RLF by Horowitz [114] and Murdock 
[115]. The nucleon form- factor model was changed be- 
tween GK by Gari and Krumpclmann [116] and the dipolc 
model. Further, the QMC model of Lu et al. [117, 118] 
predicts a density dependence for form factors that was 
calculated and applied to the GK form factors using the 
LDA (see Ref. [80]). 

Note that the calculations for the lp 3 / 2 -state in- 
clude the incoherent contributions of the unresolved 
2si/2l<^5/2-doublet. The bound-nucleon wave functions 
for these positive-parity states were taken from the 
paramctrization of Leuschner et al. and normalization 
factors were fit to said data using RDWIA calculations. 
Factors for both states of 0.12(3) relative to full occu- 
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FIG. 11: Momentum distributions for the HS, nlsh, and nlsh- 
p models. There is only a slight difference between HS and 
NLSH - for the Lp3/2-state, HS is broader spatially and thus 
drops off faster with increasing p m i BB . On the other hand, 
nlsh-p differs appreciably from both HS and nlsh, and is 
clearly distinguishable for p m i ss > 250 MeV/c for both the 
Lpi/2- and lp 3 / 2 -states. Note that both the nlsh and nlsh-p 
wave functions predict binding energies, single-particle ener- 
gies, and a charge radius for 16 O which are all in good agree- 
ment with the data. 



pancy were determined. The sensitivity of the present 
data to this incoherent admixture was evaluated by scal- 
ing the fitted doublet contribution using factors of 0.0, 
0.5, and 1.0. 

Qualitatively, the fully relativistic approach clearly did 
the best job of reproducing the data. Fully relativistic re- 
sults were shown to be much less gauge-dependent than 
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the nonrclativistic results. The CC2 current operator was 
in general less sensitive to choice of gauge, and the data 
discouraged the choice of the Weyl gauge. The differ- 
ent optical models had little effect on the shape of the 
calculations, but instead changed the overall magnitude. 
Both the GK and dipole nucleon form-factor models pro- 
duced nearly identical results. The change in the cal- 
culated GK+QMC cross section was modest, being most 
pronounced in Alt for p m iss > 300 McV/c. The results 
were best for a 100% contribution of the strength of the 
2s 1 / 2 ld 5 / 2 -doublet to the lp 3 / 2 -state, although the data 
were not terribly sensitive to this degree of freedom. 

Fig. 12 shows the left-right asymmetry Alt together 
with rdwia calculations for the removal of protons from 
the lp-shell of 16 O as a function of p m i ss for -Eboam = 
2.442 GeV. The origin of the large change in the slope of 
Alt at p m i ss « 300 MeV/c is addressed by the various 
calculations. This 'ripple' effect is due to the distortion 
of the bound-nucleon and ejectile spinors, as evidenced 
by the other three curves shown, in which the full RD- 
wia calculations have been decomposed. It is important 
to note that these three curves all retain the same basic 
ingredients, particularly the fully relativistic current op- 
erator and the upper components of the Dirac spinors. 
Of the three curves, the dotted line resulted from a cal- 
culation where only the bound-nucleon spinor distortion 
was included, the dashed line resulted from a calcula- 
tion where only the scattered-state spinor distortion was 
included, and the dashed-dotted line resulted from a cal- 
culation where undistortcd spinors (essentially identical 
to a factorized calculation) were considered. Clearly, the 
inclusion of the bound-nucleon spinor distortion is more 
important than the inclusion of the scattered-state spinor 
distortion, but both are necessary to describe the data. 

The effects of variations in the ingredients to the cal- 
culations of the left-right asymmetry Alt for the lpi/2- 
state only are shown in Fig. 13. Note that the data 
are identical to those presented in Fig. 12, as are the 
solid curves. In the top panel, the EDAI-O optical po- 
tential and nlsh bound-nucleon wave function were used 
for all the calculations, but the choice of current operator 
was varied between CCl (dashed), CC2 (solid), and CC3 
(dashed-dotted), resulting in a change in both the height 
and the p m j SS -location of the ripple in Alt- In the middle 
panel, the current operator CC2 and EDAI-O optical po- 
tential were used for all the calculations, but the choice of 
bound-nucleon wave function was varied between nlsh-p 
(dashed), nlsh (solid), and HS (dashed-dotted), result- 
ing in a change in the p m i ss -location of the ripple, but a 
relatively constant height. In the bottom panel, the cur- 
rent operator CC2 and nlsh bound-nucleon wave func- 
tion were used for all the calculations, but the choice 
of optical potential was varied between edadI (dashed), 
edai-O (solid), and edad2 (dashed-dotted), resulting in 
a change in the height of the ripple, but a relatively con- 
stant Pmisg-location. More high-precision data, particu- 
larly for 150 < p m iss < 400 MeV/c, are clearly needed 
to accurately and simultaneously determine the current 
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FIG. 12: Left-right asymmetry Alt together with rdwia cal- 
culations for the removal of protons from the lp-shell of 16 O 
as a function of p m iss for Bbeim = 2.442 GeV. Uncertainties 
are statistical (see Table XXV for the associated systematic 
uncertainties). Note that the solid curves shown here are 
identical to those shown in Figs. 13 and 15. 



operator, the bound-state wave function, the optical po- 
tential, and of course the normalization factors. This 
experiment has recently been performed in Hall A at 
Jefferson Lab by Saha et al. [119], and the results are 
currently under analysis. 



2. Comparison to rdwia, romea, and RMSGA calculations 
considering single-nucleon currents 

In this Section, the data are compared to RDWIA and 
bare romea and RMSGA calculations (which take into 
consideration single-nucleon currents only - no MEC or 
IC). The basic options employed in the calculations are 
summarized in Table VIII. Note that both the EA-bascd 
calculations stop at p m i SS = 350 McV/c as the approxi- 
mation becomes invalid. 

Fig. 14 shows measured cross-section data for the re- 
moval of protons from the lp-shell of 16 O as a function of 
p m i ss as compared to relativistic calculations at -Ebcam = 
2.442 GeV. The solid line is the RDWIA calculation, while 
the dashed and dashed-dotted lines are respectively the 
bare romea and RMSGA calculations. The normaliza- 
tion factors for the rdwia calculations are 0.73 and 0.72 
for the lpi/2-state and lp3/2-state, respectively. For the 
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FIG. 13: Left-right asymmetry Alt together with rdwia cal- 
culations for the removal of protons from the lpiy 2 -state of 
16 as a function of p m i ss for -Bbeam = 2.442 GeV. Uncertain- 
ties are statistical (see Table XXV for the associated system- 
atic uncertainties). The solid curves in all three panels are 
the same and are identical to those shown for the removal of 
protons from the lpi/2-state of le O in Figs. 12 and 15. 
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FIG. 14: Measured cross-section data for the removal of pro- 
tons from the lp-shell of 16 O as a function of p m i ss as com- 
pared to relativistic calculations at -Ebeam = 2.442 GeV. Un- 
certainties are statistical and, on average, there is an addi- 
tional ±5.6% systematic uncertainty (see Table XII) associ- 
ated with the data. The solid line is the rdwia calculation, 
while the dashed and dashed-dotted lines are respectively the 
bare ROMEA and RMSGA calculations. 



TABLE VIII: A summary of the basic options which served 
as input to the single-nucleon current rdwia, romea, and 
RMSGA comparison calculations. Results are shown in Figs. 
14 - 16. 



Input Parameter 


RDWIA 


ROMEA 
& RMSGA 


bound-nucleon wave function 


NLSH 


HS 


Optical Model 


EDAI-O 


EDAI-O 


nucleon spinor distortion 


relativistic 


relativistic 


electron distortion 


yes 


yes 


current operator 


CC2 


CC2 


nucleon form factors 


GK 


dipole 


gauge 


Coulomb 


Coulomb 



ROMEA and RMSGA calculations, they are 0.6 and 0.7 for 
the lpi/2-state and lp3/2-state, respectively. The rdwia 
calculations do a far better job of representing the data 
over the entire p m i SS range. 

Fig. 15 shows the left-right asymmetry Alt together 
with relativistic calculations for the removal of protons 
from the lp-shell of 16 O as a function of p m i ss for -Ebcam = 
2.442 GeV. The solid line is the RDWIA calculation, while 
the dashed and dashed-dotted lines are respectively the 
bare ROMEA and RMSGA calculations. Note again the 
large change in the slope of Alt at p m i ss ~ 300 MeV/c. 
While all three calculations undergo a similar change in 
slope, the rdwia calculation does the best job of repro- 



ducing the data. The romea calculation reproduces the 
data well for p m i ss < 300 MeV/c, but substantially over- 
estimates Alt for p m i ss > 300 MeV/c. The RMSGA cal- 
culation does well with the overall trend in the data, but 
struggles with reproducing the data for the T/^/2-state. 

Fig. 16 shows the Rl+tt, Rlt, and Rt effective re- 
sponse functions together with relativistic calculations 
for the removal of protons from the lp-shell of 16 as a 
function of p m iss- Note that the data point located at 
Pmiss ~ 52 MeV/c comes from the parallel kinematics 
measurements [146] (see Table XXIII), while the other 
data points come from the quasiperpcndicular kinemat- 
ics measurements (see Tables XXIV and XXV). The 
solid line is the rdwia calculation, while the dashed 
and dashed-dotted lines are respectively the bare ROMEA 
and RMSGA calculations. The agreement, particularly 
between the rdwia calculations and the data, is very 
good. The spinor distortions in the rdwia calculations 
which were required to predict the change in slope of 
Alt at p m i ss « 300 MeV/c in Fig. 12 are also essential 
to the description of Rlt- The agreement between the 
RMSGA calculations and the data, particularly for Rlt, 
is markedly poorer. 

Qualitatively, it should again be noted that none of 
the calculations presented so far have included contribu- 
tions from two-body currents. The good agreement be- 
tween the calculations and the data indicates that these 
currents are already small at Q 2 « 0.8 (GcV/c) 2 . This 
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FIG. 15: Left-right asymmetry Alt together with relativistic 
calculations of the Alt asymmetry for the removal of pro- 
tons from the lp-shell of 16 O as a function of jj m i Ha for i?beam 
= 2.442 GeV. Uncertainties are statistical (see Table XXV 
for the associated systematic uncertainties). The solid line is 
the RDWIA calculation, while the dashed and dashed-dotted 
lines are respectively the bare ROMEA and RMSGA calculations. 
Note that the solid curves shown here are identical to those 
shown in Figs. 12 and 13. 



observation is supported by independent calculations by 
Amaro et at [120, 121] which estimate the importance 
of such currents (which are highly dependent on p m i ss ) 
to be large at lower Q 2 , but only 2% for the lp^-state 
and 8% for the lp 3 / 2 -state in these kinematics. It should 
also be noted that the RDWIA results presented here are 
comparable with those obtained in independent rdwia 
analyses of our data by the Pavia Group - see Meucci et 
al. [81] 



3. Comparison to ROMEA and RMSGA calculations 
including two-body currents 

In this Section, two-body current contributions to the 
ROMEA and RMSGA calculations stemming from MEC 
and IC are presented. These contributions to the tran- 
sition matrix elements were determined within the non- 
relativistic framework outlined by the Ghent Group in 
[122, 123]. Recall that the basic options employed in the 
calculations have been summarized in Table VIII. Note 
again that both the EA-based calculations stop at p m i ss 
= 350 McV/c as the approximation becomes invalid. 

Fig. 17 shows measured cross-section data for the re- 
moval of protons from the lp-shell of 16 O as a function 
of Pmiss as compared to calculations by the Ghent Group 
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FIG. 16: Data from this work together with relativistic calcu- 
lations for the Rl+tt, Rlt, and Rt effective response func- 
tions for the removal of protons from the lp-shell of ie O as 
a function of p m iss. Uncertainties are statistical (see Tables 
XXIV, XXV, and XXIII for the associated systematic un- 
certainties). The solid line is the rdwia calculation, while 
the dashed and dashed-dotted lines are respectively the bare 
ROMEA and RMSGA calculations. 



which include MEC and IC at E hcam = 2.442 GeV. In the 
top panel, romea calculations are shown. The dashed 
line is the bare calculation, the dashed-dotted line in- 
cludes MEC, and the solid line includes both MEC and 
IC. In the bottom panel, RMSGA calculations are shown. 
The dashed line is the bare calculation, the dashed-dotted 
line includes MEC, and the solid line includes both MEC 
and IC. Note that the curves labelled 'bare' in this figure 
are identical to those shown in Fig. 14. The normal- 
ization factors are 0.6 and 0.7 for the lpi/ 2 -state and 
l.P3/2-state, respectively. The impact of the two-body 
currents on the computed differential cross section for the 
knockout of lp-shell protons from 16 O is no more than 
a few percent for low p m i ss , but gradually increases with 
increasing p m i ss - Surprisingly, explicit inclusion of the 
two-body current contributions to the transition matrix 
elements does not markedly improve the overall agree- 
ment between the calculations and the data. 

Fig. 18 shows the left-right asymmetry Alt together 
with calculations by the Ghent Group for the removal of 
protons from the lp-shell of 16 O as a function of p m i ss 
for Ebeam = 2.442 GeV. In the top two panels, romea 
calculations are shown. The dashed lines are the bare 
calculations identical to those previously shown in Fig. 
15, the dashed-dotted line includes MEC, and the solid 
line includes both MEC and IC. In the bottom panel, 
RMSGA calculations are shown. The dashed line is the 
bare calculation, the dashed-dotted line includes MEC, 
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FIG. 17: Measured cross-section data for the removal of pro- 
tons from the lp-shell of 16 O as a function of p m i ss together 
with calculations by the Ghent Group at -Ebcam = 2.442 GeV. 
Uncertainties are statistical and, on average, there is an ad- 
ditional ±5.6% systematic uncertainty (see Table XII). The 
curves labelled 'bare' are identical to those shown in Fig. 14. 



and the solid line includes both MEC and IC. While all 
three calculations undergo a change in slope at p m i ss ~ 
300 McV/c, it is again clearly the bare calculations which 
best represent the data. Note that in general, the IC were 
observed to produce larger effects than the MEC. 

Figs. 19 and 20 show the effective Rl+tt, Rlt, and 
Rt response functions together with romea and RMSGA 
calculations by the Ghent Group for the removal of pro- 
tons from the lp-shell of 16 O as a function of p m i SS - The 
dashed lines are the bare ROMEA and RMSGA calculations 
identical to those previously shown in Fig. 16, while the 
solid lines include both MEC and IC. In contrast to the 
cross-section (recall Fig. 17) and Alt (recall Fig. 18) 
situations, the agreement between the effective response- 
function data and the calculations improves with the ex- 
plicit inclusion of the two-body current contributions to 
the transition matrix elements. 



B. Higher missing energies 

In this Section, romea calculations are compared to 
the higher- E m i ss data. The basic options employed in the 
calculations have been summarized in Table VIII. 

Fig. 21 presents averaged measured cross-section data 
as a function of E m [ ss obtained at -Ebcam = 2.442 GeV for 
four discrete HRS/j angular settings ranging from 2.5° < 
8 pq < 20° , corresponding to average values of p m i ss in- 
creasing from 50 to 340 MeV/c. The cross-section values 
shown are the averaged values of the cross section mea- 
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FIG. 18: Left-right asymmetry Alt together with calcula- 
tions by the Ghent Group of the Alt asymmetry for the 
removal of protons from the lp-shell of 16 O as a function of 
Pmiss for -Ebcam = 2.442 GeV. Error bars are statistical (see 
Table XXV for the associated systematic uncertainties). The 
curves labelled 'bare' are identical to those shown in Fig. 15. 
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FIG. 19: Data from this work together with romea calcu- 
lations by the Ghent Group for the Rl+tt, Rlt, and Rt 
effective response functions for the removal of protons from 
the lp-shell of 16 O as a function of p m i ss . Uncertainties are 
statistical (see Tables XXIV, XXV, and XXIII for the associ- 
ated systematic uncertainties) . The curves labelled 'bare' are 
identical to those shown in Fig. 16. 
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FIG. 20: Data from this work together with RMSGA calcu- 
lations by the Ghent Group for the Rl+tt, Rlt, and Rt 
effective response functions for the removal of protons from 
the lp-shell of 16 O as a function of p m j aa . Uncertainties are 
statistical (see Tables XXIV, XXV, and XXIII for the associ- 
ated systematic uncertainties). The curves labelled 'bare' are 
identical to those shown in Fig. 16. 



sured on either side of q at each 9 pq . The strong peaks 
at E m i ss = 12.1 and 18.3 McV correspond to lp-shcll 
proton removal from 16 0. As in Section V A, the dashed 
curves corresponding to these peaks are the bare romea 
calculations, while the solid lines include both MEC and 
IC. The normalization factors remain 0.6 and 0.7 for the 
lpi/2- and lp3/2-states, respectively. 

For 20 < -Emiss < 30 MeV, the spectra behave in a 
completely different fashion. Appreciable strength exists 
which scales roughly with the lp-shell fragments and is 
not addressed by the present calculations of two-nucleon 
knockout. The high- resolution experiment of Leuschner 
et al. identified two additional IP3/2 fragments and sev- 
eral positive-parity states in this region which are pop- 
ulated primarily by single-proton knockout from 2p2h 
components of the ground-state wave function. Two- 
body currents and channel-coupling in the final state 
also contribute. This strength has also been studied 
in (7,p) experiments, and has been interpreted by the 
Ghent Group [124] as the post-photoabsorption popu- 
lation of states with a predominant \p2h character via 
two-body currents. 

For -E m i SS > 30 MeV, in the top panel for p m iss = 50 
MeV/c, there is a broad and prominent peak centered at 
E m iss ~ 40 MeV corresponding largely to the knockout of 
ls^-state protons. As can be seen in the lower panels, 
the strength of this peak diminishes with increasing p m iss , 
and completely vanishes beneath a flat background by 
p miss = 280 McV/c. For E miss > 60 MeV and p miss > 



280 MeV/c, the cross section decreases only very weakly 
as a function of p m iss, and is completely independent of 
E ■ 

-^miss • 

In order to estimate the amount of the cross section 
observed for Emiss > 25 MeV that can be explained by 
the single-particle knockout of protons from the lsi/2- 
state, the data were compared to the romea calcula- 
tions of the Ghent Group. The dashed curves are the 
bare calculations, while the solid lines include both MEC 
and IC. A normalization factor of 1.0 for the lsj/2-state 
single-particle strength was used. The two calculations 
are indistinguishable for p m iss < 145 MeV/c, and the 
agreement between these calculations and the measured 
cross-section data is reasonable (see the top two panels of 
Fig. 21 where there is an identifiable ls^-state peak at 
E m iss ~ 40 MeV) . At higher p m ; ss (where there is no clear 
lsx/2-state peak at E m iss ~ 40 MeV), the data are sub- 
stantially larger than the calculated bare cross section. 
Inclusion of MEC and IC improves the agreement, but 
there is still roughly an order-of-magnitudc discrepancy. 
The RDWIA calculations demonstrate similar behavior. 
Thus, the p m i ss > 280 MeV/c data are not dominated 
by single-particle knockout. Note that the magnitude of 
(St — Sl) is consistent with that anticipated based on the 
measurements of Ulmer et al. at Q 2 = 0.14 (GeV/c) 2 and 
Dutta et al. at Q 2 = 0.6 and 1.8 (GcV/c) 2 . Together, 
these data suggest that transverse processes associated 
with the knockout of more than one nucleon decrease 
with increasing Q 2 . 

Also shown as dashed-dotted curves in Fig. 21 are the 
calculations by the Ghent Group [125] for the (e, e'pp) 
and (e, e'pn) contributions to the (e, e'p) cross section 
performed within a Hartree-Fock framework. This two- 
particle knockout cross section was determined using 
the Spectator Approximation, in a calculation which in- 
cluded MEC, IC, and both central short-range correla- 
tions (SRC) and tensor medium-range correlations. Note 
that in these kinematics, this calculation performed with 
SRC alone produced only 2% of the two-particle knock- 
out cross section, while including both SRC and ten- 
sor correlations produced only 15% of the two-particle 
knockout cross section. The calculated two-particle 
knockout cross section is essentially transverse in nature, 
since the two-body currents are predominantly trans- 
verse. The calculated strength underestimates the mea- 
sured cross section by about 50% but has the observed 
flat shape for E m iss > 50 MeV. It is thus possible that 
heavier meson exchange and processes involving three (or 
more) nucleons could provide a complete description of 
the data. 

The measured effective response functions R l+ tt, 
Rlt, and Rt together with romea calculations for p m ; ss 
= 145 MeV/c and p m iss = 280 MeV/c are presented in 
Fig. 22. Kinematic overlap restricted separations to 
Emiss < 60 MeV. The dashed curves are the bare romea 
calculations, while the solid curves include both MEC 
and IC. Also shown as dashed-dotted curves are the inco- 
herent sum of these 'full' calculations and the computed 
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FIG. 21: Data from this work together with ROMEA calculations by the Ghent Group for the i? m i ss -dependence of the cross 
section obtained at -Bbeam = 2.442 GeV. The data are the averaged cross section measured on either side of q at each 9 pq . 
Normalization factors of 0.6, 0.7, and 1.0 have been used for the lpi/2-, lf>3/2-> an d lsi/2-states, respectively. Uncertainties are 
statistical and, on average, there is an additional ±5.9% systematic uncertainty (see Tables XVIII - XXII) associated with the 
data. Also shown are calculations by the Ghent Group for the (e, e'pN) contribution. 



(e, e'pN) contribution. In general, the data do not show 
the broad peak centered at -E m i ss ~ 40 MeV correspond- 
ing to the knockout of ls^-state protons predicted by 
the calculations. At p m iss = 145 MeV/c, the bare cal- 
culation is consistently about 60% of the magnitude of 
the data. Inclusion of MEC and IC does not apprecia- 
bly change the calculated Rl+tt, but does improve the 
agreement between data and calculation for Rlt and Rt- 



The measured response Rl+tt (which is essentially equal 
to R L since ^-Rtt is roughly 7% of Rl in these kine- 
matics - see Ref. [80]) is larger than the calculation for 
-Emiss < 50 MeV and smaller than the calculation for 
Smiss > 50 MeV. The agreement between the calculation 
and the data for Rlt is very good over the entire E m [ SS 
range. Since the measured response Rlt is nonzero for 
^miss > 50 MeV, the measured response Rl must also be 
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FIG. 22: Data from this work together with ROMEA calcu- 
lations by the Ghent Group for the _E m i aa -dependence of the 
Rl+tt, Rlt, and Rt effective response functions. Uncer- 
tainites are statistical (see Tables XXVII, XXVIII, and XXVI 
for the associated systematic uncertainties). Also shown is the 
(e, e'pN) contribution. 



nonzero. The measured response Rt is somewhat larger 
than the calculation for i? m i ss < 60 MeV. 

At pmiss = 280 MeV/c, the bare calculation does not 
reproduce the .Emiss-depcndence of any of the measured 
effective response functions. The inclusion of MEC and 
IC in the calculation substantially increases the magni- 
tude of all three calculated response functions, and thus 
improves the agreement between data and calculation. 
The measured Rl+tt (which is dominated by Rl) is 
consistent with both the calculation and with zero. The 
measured Rlt is about twice the magnitude of the cal- 
culation. Since the measured Rlt is nonzero over the en- 
tire E m i ss range, the measured Rl must also be nonzero. 
The measured Rt is significantly larger than both the 
calculations and nonzero out to at least -E m iss ~ 60 MeV. 
The fact that the measured Rt is much larger than the 
measured Rl indicates the cross section is largely due to 
transverse two-body currents. And finally, it is clear that 
(e, e'pN) accounts for a fraction of the measured trans- 
verse strength which increases dramatically with increas- 
ing p miss . 

Fig. 23 shows the calculations by the Ghent Group 
[105] of the contribution to the differential 16 0(e,e'p) 
cross section from two-nucleon knockout as a function 
of E miss and 9 p for -Ebcam = 2.442 GeV. The upper- 
left panel shows the contribution of central correlations. 
The upper-right panel shows the combined contribution 
of central and tensor correlations. Tensor correlations 
are anticipated to dominate central correlations over the 



ranges of E m iss and p m iss investigated in this work. The 
lower-left panel shows the combined contribution of cen- 
tral and tensor correlations (two-nucleon correlations) to- 
gether with MEC and IC (two-body currents). Two-body 
currents are anticipated to dominate two-nucleon corre- 
lations over the ranges of -E m iss and p m ; ss investigated in 
this work. For convenience, the variation of p m i ss with 
.Emiss and 9 P is shown in the bottom-right panel. 



VI. KINEMATIC CONSISTENCY OF lp-SHELL 
NORMALIZATION FACTORS 



There has been longstanding discussion regarding the 
reliability of the spectroscopic factors determined for dis- 
crete states from single- nucleon electromagnetic knock- 
out. Recently, there has been speculation that these fac- 
tors might appear to increase with Q 2 as a quasiparticle 
state is probed with finer resolution. In this Section, 
Kelly [126] has used the RDWIA to analyze the normal- 
ization factors fitted to the available 16 0(e, e'p) data for 
the lpi/2- and lp 3 / 2 -states obtained in the experiments 
summarized in Table IX. If the rdwia model is accurate, 
these factors should be independent of the experimental 
kinematics. 



A. RDWIA analysis of the available le O(e, e'p) 
data 

The rdwia calculations used in this procedure cm- 
ployed the f 2 off-shell single-nucleon current operator 
with the MMD form factors of Mergell et al. [127] in 
the Coulomb gauge. The partial- wave expansions were 
performed using the second-order Dirac equation, includ- 
ing spinor distortion. Nucleon distortion was evaluated 
using the EDAI-O and edadI optical potentials, and elec- 
tron distortion was evaluated in the (? e ff Approximation. 
No attempt to directly fit the overlap functions to the 
knockout data has been made here. Instead, the recently 
developed wave functions HS, nlsh, and nlsh-p (recall 
Fig. 11) were again tested, this time to see if they could 
satisfactorily reproduce the experimental p m i ss distribu- 
tions independently of Q 2 . 

The results are expressed in terms of normalization fac- 
tors which compare a RDWIA calculation for a fully occu- 
pied subshell with experimental data and are presented 
in Table X. These factors were obtained by least-squares 
fitting to the data in the range p m i SS < 200 MeV/c where 
the RDWIA should be most reliable. When experimen- 
tally unresolved, the contamination of the lp 3 / 2 -state by 
the ld5/22s 1 / 2 -doublet was included by the incoherent 
summation of the paramctrizations of Lcuschncr et al. 
as previously described. 

The data sets demonstrated a slight preference for the 
EDADl optical potential over the EDAI-O optical poten- 
tial. This was concluded based on the quality of the fits 
and the more consistent nature of the extracted normal- 
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FIG. 23: (Color online) Calculations by the Ghent Group of the contribution to the differential 1 0(e, e'p) cross section from 
two-nucleon knockout as a function of -E m i ss and 9 P for Bbeam = 2.442 GeV. The upper-left panel shows the contribution of 
central correlations. The upper-right panel shows the combined contribution of central and tensor correlations. The lower-left 
panel shows the combined contribution of central and tensor correlations (two-nucleon correlations) together with MEC and 
IC (two-body currents). The relationship between the various kinematic quantities is shown in the bottom-right panel. 



TABLE IX: A summary of the kinematic conditions for the data examined in the 16 0(e,e'p) consistency study. 



label 


authors 


kinematics 


T P 
(MeV) 


Q 2 

(GeV/c) 2 


X 


2s!/ 2 la!5/2-doublet 


data 


a 
b 
c 


Leuschner et al. [14] 
Spaltro et al. [13] 
Chinitz et al. [12] 


parallel 
perpendicular 
perpendicular 


96 
84 
160 


varied 
0.20 
0.30 


varied 
1.07 
0.91 


resolved 
resolved 
computed a 


reduced a 
differential a 
differential cr 


d 


this work 


perpendicular 


427 


0.80 


0.96 


computed 6 


differential cr 


e 
f 

g 


Bernheim et al. [11] 
Blomqvistl et al. [15] 
Blomqvist2 et al. [15] 


perpendicular 
parallel 
highly varied 


100 
92 
215 


0.19 
0.08 
0.04 - 0.26 


0.90 
0.30 - 0.50 
0.07 - 0.70 


computed b 
resolved 
resolved 


reduced cr 
reduced cr 
reduced cr 



"The lp 3 / 2 -state data were corrected for the contamination of the 1^5/2 2si/2 _ doublet by Chinitz et al. 

b The contamination of the lp3/2-state by the ld5/22si/ 2 -doublet was computed according to the method outlined in Section V A 1. 



ization factors for low Q 2 . None of the variations consid- 
ered in Table X (nor any of those considered in Table VII 
for that matter) were able to reproduce the bp 3 / 2 -state 
for data set (b) in the range 50 < p m iss < 120 MeV/c. 
This problem is also responsible for the discrepancy seen 
in Fig. 1 for R LT at Q 2 = 0.2 (GcV/c) 2 , and has not yet 



been explained satisfactorily. 

Unfortunately, none of the selected wave functions pro- 
vided an optimal description of the experimental p m ; ss 
distributions independent of Q 2 . Fig. 24 shows a sam- 
ple set of fits to the various 16 0(e,e'p) data sets based 
on the HS bound- nucleon wave function and the edadI 
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TABLE X: Normalization factors deduced for the data sets presented in Table IX for p m i ss < 200 MeV/c. The first term in 
each column is for the Lpi/2-state, while the second term is for the lp3/2-state. 





EDAI-O 


EDADl 


HS 

s a x 2 


NLSH 

s a x 2 


NLSH-P 

s a x 2 


HS 

s a x 2 


NLSH 

s a x 2 


NLSH-P 

s a x 2 


a 
b 

c 


0.55 0.46 
0.61 0.66 
0.54 0.56 


0.9 4.2 
2.7 6.3 
8.7 17.9 


0.60 0.47 
0.65 0.68 
0.60 0.58 


2.5 6.0 
5.4 8.0 
24.8 25.3 


0.53 0.41 
0.58 0.58 
0.51 0.47 


1.0 1.5 
2.3 3.6 
9.0 15.0 


0.60 0.55 
0.71 0.75 
0.59 0.61 


0.8 2.4 
2.3 5.1 
8.0 18.2 


0.66 0.56 
0.77 0.77 
0.66 0.63 


2.3 3.7 
4.2 6.4 
16.7 22.1 


0.57 0.48 
0.68 0.65 
0.56 0.50 


0.8 1.4 
2.2 4.1 
7.2 23.8 


(1 


0.62 0.63 


30.8 4.7 


0.70 0.65 


0.5 6.7 


0.59 0.52 


19.6 15.3 


0.62 0.63 


32.5 2.4 


0.70 0.66 


1.2 3.9 


0.60 0.53 


20.0 14.4 


e 
f 
g 


0.43 0.46 
0.53 0.41 
0.42 0.37 


1.0 1.9 
3.0 4.2 
2.0 1.4 


0.48 0.47 
0.54 0.42 
0.44 0.37 


2.2 2.4 
5.0 5.8 
4.7 2.0 


0.42 0.40 
0.51 0.38 
0.40 0.33 


1.0 1.1 
2.7 1.9 
2.5 5.9 


0.48 0.52 
0.57 0.50 
0.42 0.40 


1.0 1.5 
3.2 3.7 
1.8 1.9 


0.54 0.54 
0.59 0.51 
0.44 0.41 


1.5 1.9 
5.8 5.0 

6.6 2.8 


0.47 0.45 
0.54 0.46 
0.40 0.36 


1.0 1.3 
2.7 2.1 
1.9 5.4 



optical potential. From a x 2 perspective, it is clear that 
NLSH offered the best description of the data at Q 2 = 0.8 
(GeV/c) 2 , but that most of the lower Q 2 data are best de- 
scribed by either NLSH-P or nlsh; in fact, HS may be the 
best compromise currently available. Data set (g) from 
Mainz suggested a substantially different normalization. 

An estimate of the uncertainty in the normalization 
factors due to variations of the bound- nucleon wave func- 
tion was made by comparing nlsh/hs and nlsh-p/hs 
normalization-factor ratios for each state holding the op- 
tical potential and other model input constant. For the 
lowest-lying lpi/2- and lp3/ 2 -states and averaged over 
all data sets, nlsh/hs ratios of 1.08 and 1.02 were ob- 
tained. Similarly, nlsh-p/hs ratios of 0.96 and 0.87 
were obtained. These ratios are qualitatively consis- 
tent with the behavior of the data and the calculations 
near the peaks of the momentum distributions shown in 
Fig. 11. Therefore, a cautious estimate of the uncer- 
tainty due to the bound-nucleon wave function is of order 
±10%. Further, by changing the optical potentials be- 
tween edai-O and EDADl and holding the bound-nucleon 
wave function and other model input constant, edai- 
o/edadI normalization-factor ratios for a given data set 
with Q 2 < 0.4 (GeV/c) 2 averaged to about 0.90. This ra- 
tio became 0.98 at Q 2 = 0.8 (GeV/c) 2 , where the attenu- 
ation in the potentials is practically identical. Therefore, 
the uncertainty due to variations of the optical poten- 
tial is at least ±5% and would probably be larger if the 
sample of 'reasonable' potentials were expanded. 

The information presented in Table VII suggests that 
there would be similar uncertainties in the normalization 
factors for the Q 2 = 0.8 (GeV/c) 2 data arising from Gor- 
don and gauge ambiguities in the single-nucleon current 
operator. Note that values of S a for the same model are 
generally larger in Table VII than in Table X because 
the former summarizes a study of the entire p m i ss range 
while the latter is limited to p m i SS < 200 MeV/c, where 
the reaction model is likely to be most accurate. Data 
for larger p m iss also tend to have a higher % 2 . These 
problems for large p m iss may arise from inaccuracies in 
the bound-nucleon wave functions above the Fermi mo- 
mentum, neglecting two-body currents, neglecting chan- 
nel coupling in the final state, or density-dependence in 



the form factors, to name a few. Therefore, a realis- 
tic estimate of the model dependence of the normaliza- 
tion factors for (e,e'p) reactions should not be less than 
±15%. This estimated precision is consistent with that 
suggested in Rcf. [2] - although the relativistic model 
improves our description of Alt, recoil polarization, and 
other normalization-independent features of the reaction, 
the model dependencies that affect the normalization un- 
certainty are not significantly improved. Assuming that 
the reaction model is most reliable at large Q 2 and mod- 
est Pmiss, the six lpi/2- and Lp 3 / 2 -statc normalization 
factors for data set (d) in Table X (this work) were av- 
eraged to conclude that the normalization factors for the 
lowest lpi/2- and lp 3 /2-states in 15 N are approximately 
0.63(9) and 0.60(9). 

As previously mentioned, Leuschncr et al. identified 
two additional lp 3 / 2 -states with excitation energies be- 
tween 9 and 13 MeV that together carry approximately 
11% of the strength of the lowest-energy fragment. How- 
ever, those states were not resolved by the present ex- 
periment. If the assumption is made that the same ra- 
tio applies at Q 2 ~ 0.8 (GeV/c) 2 , then the total I-P3/2 
strength below 15 MeV excitation is estimated to be ap- 
proximately 67% of full occupancy, and the total lp-shell 
spectroscopic strength below 15 MeV represents 3.9 ±0.6 
protons or about 65% of full occupancy. This result re- 
mains 10—20% below predictions from recent calculations 
of the hole spectral function by Barbieri and Dickhoff 
[128], but no experimental estimate for the additional p- 
shell strength that might lurk beneath the continuum is 
available. 

To obtain more precise normalization factors, it would 
be necessary to apply a relativistic analysis to data in 
quasipcrpcndicular kinematics for several values of Q 2 
larger than about 0.5 (GeV/c) 2 and with sufficient cov- 
erage of the Pmiss distribution to fit the wave function, 
requiring that the bound-nucleon wave function be in- 
dependent of Q 2 . Although such data do not yet exist 
for 16 O, the recently completed experiment in Hall A at 
Jefferson Lab will provide substantially more data points 
for the critical p m iss < 200 MeV/c region for Q 2 = 0.9 
(GeV/c) 2 . 
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FIG. 24: Fits to various le O(e, e'p) data sets based on the HS bound-nucleon wave function and the edadI optical potential. See 
Table IX for the key to the data-set labels. Open points and solid lines pertain to the lp 1( / 2 -state, while solid points and dashed 
lines pertain to the lp 3 / 2 -state. The dashed-dotted lines include the contributions of the positive parity 2si/2lds/2-doublet to 
the lp 3 / 2 -state. Panel (d) shows the data from this work. 



B. Q 2 -dependence of normalization factors 

Lapikas et al. [129] have performed a similar type of 
analysis of the (5 2 -dcpendence of the normalization fac- 



tors for the 12 C(e, e'p) reaction. In their work, several 
data sets with Q 2 < 0.3 (GeV/c) 2 were analyzed using 
a nonrelativistic dwia model. For each data set, a nor- 
malization factor and the radius parameter for a Woods- 



28 



Saxon binding potential were fitted to the reduced cross 
section for discrete states, and the potential depths were 
adjusted to fit the separation energies. Consistent nor- 
malization factors were obtained for all data sets save 
those measured by Blomqvist et al. at Mainz [130]. A 
new experiment was thus performed at NIKHEF du- 
plicating the Mainz kinematics. The new results were 
also consistent with all data sets save those from Mainz. 
Lapikas et al. thus concluded that the Mainz data were 
normalized incorrectly. Recall that similar doubts re- 
garding the normalization of the companion 16 0(e,e'p) 
experiment [15] at Mainz had been expressed earlier by 
Kelly [4] , but independent data duplicating the measure- 
ment are unfortunately not available. After excluding the 
Mainz data, Lapikas et al. determined that the summed 
lp-shcll strength for 12 C could in fact be deduced from 
data for Q 2 < 0.3 (GeV/c) 2 with an uncertainty of ±3%. 
However, they did not consider the effects of variations 
of the optical model or several other uncertain aspects 
of the reaction model. As discussed previously, a more 
realistic estimate of the relative uncertainty in the nor- 
malization factors must be closer to ±15% due to the 
inevitable model dependence of the dwia. Furthermore, 
it is possible that variation of the Woods-Saxon radius 
might affect the resulting normalization factors. If the 
overlap function is an intrinsic property of the nuclear 
wave function, it should not depend upon Q 2 . Further, 
it should be possible to fit a common radius to all data 
simultaneously; if not, the accuracy of the reaction model 
must be questioned. And of course, it has been demon- 
strated in recent years that a relativistic DWIA model is 
preferable to a nonrclativistic approach. 

Lapikas et al. also used the bound-nucleon wave func- 
tions and normalization factors obtained from their non- 
rclativistic analysis at low Q 2 to analyze the transparency 
of 12 C for Q 2 up to 7 (GeV/c) 2 and the summed lp- 
and ls-shell spectroscopic amplitude. They found the 
summed spectroscopic strength was approximately con- 
stant at 0.58 for Q 2 < 0.6 (GeV/c) 2 , but rose for larger 
Q 2 and appeared to approach the Independent-Particle 
Model limit of unity somewhere near Q 2 « 10 (GeV/c) 2 . 
They speculated that the apparent Q 2 -dependence of this 
spectroscopic strength might be related to the resolution 
at which a quasiparticle is probed, with long-range cor- 
relations that deplete the single-particle strength becom- 
ing less important at higher Q 2 and finer resolution. A 
subsequent analysis by Frankfurt et al. using Glauber 
calculations for heavier targets [131] supports this inter- 
pretation. 

Little evidence is seen here for a systematic dependence 
in the normalization factors upon cither T p or Q 2 for the 
lowest lp-states of 16 O for the data that are presently 
available. Unfortunately, these data do not reach high 
enough Q 2 to address the resolution hypothesis. Further- 
more, the normalization factors for two of the data sets 
appear to be anomalously low. A normalization problem 
might not be too surprising for data set (e) because it 
comes from one of the earliest experiments on this reac- 



tion, but data set (g) comes from a fairly recent experi- 
ment at Mainz and uses an ejectile energy large enough 
for the reaction model to be reliable. As discussed above 
for the case of 12 C, it is likely that data set (g) also has 
a normalization error [147]. If these two data sets are 
disregarded, the remaining low Q 2 data are consistent 
with the normalization factors deduced from the current 
Q 2 = 0.8 (GeV/c) 2 data. 



VII. SUMMARY AND CONCLUSIONS 

The 16 0(e,e'p) reaction in QE, constant (q, ui) kine- 
matics at Q 2 « 0.8 (GeV/c) 2 , q « 1 GeV/c, and w « 445 
McV was studied for < E miss < 120 McV and < p raiss 

< 350 MeV/c. Five-fold differential cross-section data for 
the removal of protons from the lp-shcll were obtained 
for < p m i ss < 350 McV/c. Six- fold differential cross- 
section data for < E miss < 120 MeV were obtained for 

< Pmiss < 340 MeV/c. These results were used to extract 
the Alt asymmetry and the Rl, Rt, Rl+tt, and Rlt 
effective response functions over a large range of E m [ ss 
and p miss . 

The data were interpreted in subsets corresponding 
to the lp-shcll and the ls^-state and continuum, re- 
spectively, lp-shell data were interpreted within three 
fully relativistic frameworks for single-particle knockout 
which do not include any two-body currents: rdwia, 
romea, and RMSGA. Two-body current contributions 
to the ROMEA and RMSGA calculations for the lp-shell 
stemming from MEC and IC were also considered. The 
lS]/2-state and continuum data were considered within 
the identical ROMEA framework both before and after 
two-body current contributions due MEC and IC were 
included, (e, e'pN) contributions to these data were also 
examined. 

Overall, the rdwia calculations provided by far the 
best description of the lp-shell data. Dynamic effects 
due to the inclusion of the lower components of the 
Dirac spinors in these calculations were necessary to self- 
consistcntly reproduce the lp-shcll cross-section data, the 
Alt asymmetry, and the Rlt effective response function 
over the entire measured range of p m iss- Within the RD- 
WIA framework, the four most important ingredients were 
the inclusion of both bound-nucleon and ejectile spinor 
distortion, the choice of current operator, the choice of 
bound-nucleon wave function, and the choice of optical 
potential. Inclusion of the spinor distortion resulted in a 
diffractive change in slope in Alt at p m i ss ~ 300 MeV/c 
which agreed nicely with the data. A different choice 
of current operator either damped out or magnified this 
change in slope. A different choice of bound-nucleon 
wave function changed the p m is S -location of the change 
in slope, but preserved the magnitude. A different choice 
of optical potential changed the magnitude of the change 
in slope but preserved the p m i ss -location. 

As anticipated, since p p ss 1 GeV/c, the romea calcu- 
lations provided a reasonable description of the lp-shcll 
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data. For this energy range, optical models generally pro- 
vide an overall better description of proton elastic scat- 
tering than does the Glauber model. This is in part due 
to important medium modifications of the NN interac- 
tion from Pauli blocking and spinor distortion. Surpris- 
ingly, the unfactorized 'out-of-the-box' RMSGA calcula- 
tion provided a fairly good description of the lp-shcll 
data already at this relatively low proton momentum. 
Adding the contributions of two-body currents due to 
MEC and IC to the descriptions of the lp-shell data pro- 
vided by the bare romea and RMSGA calculations did 
not improve the agreement. 

The rdwia calculation with single-nucleon currents 
was used to fit normalization factors to the data from 
this experiment and from several other experiments at 
lower Q 2 . Ignoring two experiments which appear to have 
normalization problems, normalization factors of 0.63(9) 
and 0.60(9) were obtained for the lowest lpi/2- and IP3/2- 
states with no significant dependence upon Q 2 or T p . The 
estimated uncertainties account for variations due to the 
choice of bound-nucleon wave functions, optical poten- 
tials, and other aspects of the model. After accounting 
for other known but unresolved lp 3 / 2 -states, the total 
lp-shcll spectroscopic strength below about 15 MeV ex- 
citation is estimated to be about 0.65 ± 0.10 relative to 
full occupancy. 

For 25 < E miss < 50 MeV and p miss < 145 McV/c, 
the reaction was dominated by the knockout of lS]/2- 
state protons and the cross section and effective response 
functions were reasonably well-described by bare ROMEA 
calculations which did not consider the contributions of 
two-body currents due to MEC and IC. However, as p m iss 
increased beyond 145 MeV/c, the single-particle aspect of 
the reaction diminished. Cross-section data and response 
functions were no longer peaked at E m i ss ss 40 MeV, nor 
did they exhibit the Lorentzian s-shell shape. Already 
at Pmiss = 280 MeV/c, the same bare romea calcula- 
tions that did well describing the data for p m iss < 145 
MeV/c underestimated the cross-section data by more 
than an order of magnitude. Including the contributions 
of two-body currents due to MEC and IC improved the 
agreement for -E m iss < 50 MeV, but the calculations still 
dramatically undcrpredict the data. 

For 25 < E miss < 120 MeV and p miss > 280 MeV/c, 
the cross-section data were almost constant as a func- 
tion of both Pmiss and -E m iss- Here, the single-particle 
aspect of the ls^-state contributed <10% to the cross 
section. Two-nuclcon (e, e'pN) calculations accounted 
for only about 50% of the magnitude of the cross-section 
data, but reproduced the shape well. The model, which 
explained the shape, transverse nature, and 50% of the 
measured cross section, suggested that the contributions 
of the two-nucleon currents due to MEC and IC are much 
larger than those of the two-nucleon correlations. The 
magnitude of the measured cross section that remains un- 
accounted for suggests additional currents and processes 
play an equally important role. 
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APPENDIX A: QUASIELASTIC RESULTS 



1. Cross-section data 



The cross-section data are presented in this section. Note that the data sets obtained for 9 pq = 0° presented in 
Tables XI, XIII, XVI, and XVIII do not truly represent parallel kinematics even though the HRS/j was aligned along q. 
Since p p and q had about the same magnitude, p m j ss arose from the slight angles between them, not from differences 
in their magnitudes. However, since the distribution of p miss was symmetrical about q, the conditions of parallel 
kinematics were closely approximated. 



a. lp-shell 



Cross-section data for QE proton knockout from the lp-shell of 16 are presented in Tables XI and XII. 

TABLE XI: Measured cross-section data for QE proton knockout from the lp-shcll 
of le O for < T p > = 427 MeV. The p mis3 bins were 20 MeV/c wide. Cuts were 
applied to remove the radiative tail from 1 H(e, ep) such that < p m i ss > — 52.5 
McV/c in each case. 



-fib cam pq 

(Gcv) n 


lPl/2 -s t, a t. c 


lp 3( / 2 -state 


< Q A > < u > d°a/duidn c dn p (sys) 
(GcV/c) 2 (MeV) (nb/MeV/sr 2 ) (%) 


< Q A > < uj > d°a/dLodn,,dn p (sys) 
(GcV/c) 2 (MeV) (nb/MeV/sr 2 ) (%) 


0.843 0.0 


0.810 436.0 0.0922 ± 0.0118 5.4 


0.810 436.0 0.1143 ± 0.0115 5.4 


1.643 0.0 


0.820 421.5 0.5827 ± 0.0486 5.5 


0.820 421.5 0.7418 ± 0.0514 5.4 


2.442 0.0 


0.815 423.0 1.5030 ± 0.1380 5.5 


0.815 423.0 1.8540 ± 0.1500 5.5 



TABLE XII: Measured cross-section data for QE proton knockout from the lp-shcll 
of ls O for < Q 2 > = 0.800 (GeV/c) 2 , < u > = 436 MeV, and < T p > = 427 MeV. 
The Pmiss bins were 20 MeV/c wide. 









lPl/2 


state 




1P3/2 


state 


-Kb cam 




< Pmiss > 


d tr / d ui d O e d £1 v 


(sys) 


< Pmiss > 


d° a /dudQ G dn p 


(sys) 


(GcV) 


C) 


(McV/c) 


(nb/MeV/sr 2 ) 


(%) 


(McV/c) 


(nb/MeV/sr 2 ) 


(%) 


0.843 


8.0 


150.0 


0.0789 ± 0.0057 


5.6 


150.0 


0.1467 ± 0.0076 


5.5 




16.0 


275.0 


0.0011 ± 0.0003 


6.1 


275.0 


0.0058 ± 0.0006 


7.2 


1.643 


-8.0 


-148.0 


0.2950 ± 0.0320 


5.5 


-146.0 


0.5160 ± 0.0370 


5.5 




8.0 


148.0 


0.5250 ± 0.0310 


5.5 


146.0 


1.0390 ± 0.0360 


5.4 


2.442 


-20.0 


-355.0 


0.0023 ± 0.0011 


5.5 


-355.0 


0.0054 ± 0.0011 


5.4 




-16.0 


-279.0 


0.0143 ± 0.0029 


5.7 


-275.0 


0.0288 ± 0.0051 


6.1 




-8.0 


-149.0 


0.9060 ± 0.0260 


5.5 


-149.0 


1.5740 ± 0.0374 


5.5 




-2.5 


-60.0 


1.5981 ± 0.0456 


5.4 


-60.0 


2.2360 ± 0.0540 


5.5 




2.5 


60.0 


1.5380 ± 0.0513 


5.4 


60.0 


2.6210 ± 0.0650 


5.5 




8.0 


149.0 


1.4605 ± 0.0261 


5.5 


149.0 


2.9950 ± 0.0374 


5.5 




16.0 


279.0 


0.0303 ± 0.0029 


5.7 


276.0 


0.1672 ± 0.0051 


0.2 




20.0 


330.0 


0.0057 ± 0.0005 


5.6 


330.0 


0.0214 ± 0.0008 


5.5 



b. Higher missing energies 

Cross-section data together with statistical uncertainties for QE proton knockout from 16 O for .E^ss > 20 MeV are 
presented in Tables XIII - XXII. 

TABLE XIII: Measured cross-section data for QE proton knockout from ie O for 
Ebeam = 0.843 GeV, 6 pq = 0.0°, and E miss > 25 MeV. The p mies bins were 5 
McV/c wide. Cuts were applied to remove the radiative tail from 1 H(e, ep). There 
is an additional 5.8% systematic uncertainty associated with these results. 



fimisa 


< w > 


< Q 2 > 


< Pmiss > 


d° cr j ' du)dE p d£l e di~l p 


(MeV) 


(MeV) 


(GeV/c) 2 


(McV/c) 


(nb/MeV 2 /sr 2 ) 


27.5 


445.2 


0.796 


38.7 


0.0024 ± 0.0005 


32.5 


445.2 


0.795 


41.8 


0.0085 ± 0.0008 


37.5 


445.8 


0.794 


45.0 


0.0097 ± 0.0008 


42.5 


446.3 


0.793 


48.5 


0.0113 ± 0.0009 


47.5 


447.7 


0.790 


51.3 


0.0106 ± 0.0010 


52.5 


449.7 


0.786 


53.2 


0.0065 ± 0.0010 


57.5 


451.7 


0.782 


55.5 


0.0062 ± 0.0013 
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TABLE XIV: Measured cross-section data for QE proton knockout from O for 
E hcam = 0.843 GcV, 6 pq = 8.0°, and E miss > 25 MeV. The p miss bins were 5 
McV/c wide. There is an additional 6.0% systematic uncertainty associated with 
these results. 





< to > 


< 


Q 2 > 


< Pmi 


» > 




(McV) 


(McV) 


(GcV/c) 2 


(McV/c) 


(nb/McV 2 /sr 2 ) 


25.5 


444.9 





796 


143 


3 





0048 


b 


0.0013 


26.5 


444.6 





796 


143 


5 





0019 


i 


0.0010 


27.5 


445.2 





795 


143 


8 





0008 


i 


0.0009 


28.5 


444.8 





796 


143 


5 





0041 


i 


0.0011 


29.5 


444.6 





796 


143 


9 





0024 


± 


0.0010 


30.5 


444.7 





796 


143 


9 





0022 


± 


0.0009 


31.5 


444.9 





796 


144 


2 





0022 


i 


0.0009 


32.5 


444.8 





796 


144 








0041 


i 


0.0010 


33.5 


445.0 





795 


144 


6 





0023 


± 


0.0009 


34.5 


445.0 





795 


144 


5 





0013 




0.0008 


35.5 


444.9 





796 


144 


7 





0031 


± 


0.0009 


36.5 


445.0 





796 


145 








0018 


± 


0.0008 


37.5 


445.1 





796 


145 


1 





0042 


± 


0.0010 


38.5 


444.9 





796 


145 


9 





0013 


b 


0.0008 


39.5 


444.9 





796 


145 


6 





0038 


i 


0.0010 


40.5 


445.0 





796 


146 


2 





0027 


± 


0.0009 


41.5 


444.9 





796 


146 


2 





0030 


b 


0.0008 


42.5 


445.1 





795 


146 


4 





0029 


i 


0.0009 


43.5 


445.1 





795 


147 


3 





0016 


;b 


0.0008 


44.5 


444.6 





796 


147 


1 





0020 


I 


0.0008 


45.5 


445.4 





795 


147 


6 





0023 


± 


0.0008 


46.5 


446.1 





793 


147 


8 





0028 




0.0009 


47.5 


446.5 





793 


148 








0040 


1 


0.0010 


48.5 


447.1 





791 


148 


4 





0016 


b 


0.0009 


49.5 


447.5 





790 


148 


5 





0033 


'I 


0.0010 


50.5 


448.0 





789 


149 


3 





0029 


± 


0.0010 


51.5 


448.5 





789 


149 


5 





0015 




0.0009 


52.5 


449.0 





787 


149 


8 





0013 


i 


0.0007 


53.5 


449.4 





786 


150 








0024 


± 


0.0007 


54.5 


450.0 





785 


151 


1 





0022 


"I 


0.0009 


55.5 


450.5 





784 


151 


2 





0008 


i 


0.0008 


56.5 


451.1 





783 


151 


7 





0019 


i 


0.0010 


57.5 


451.6 





782 


151 


7 





0010 


1 


0.0010 


58.5 


452.0 





782 


152 


7 





0007 


1 


0.0009 


59.5 


452.5 





781 


152 


6 





0016 


1 


0.0011 



TABLE XV: Measured cross-section data for QE proton knockout from O for 
E beam = 0.843 GeV, 6 vq = 16.0°, and E miss > 25 McV. The p miss bins were 5 
McV/c wide. There is an additional 5.9% systematic uncertainty associated with 
these results. 





<u > 


< Q 2 > 


< Pmiss > 


d° a I ' dujdE p dQ. e di1 v 


(McV) 


(McV) 


(GcV/c) 2 


(McV/c) 


(nb/McV 2 /sr 2 ) 


27.5 


446.0 


0.795 


281.6 


0.0004 ± 0.0001 


32.5 


446.0 


0.795 


281.3 


0.0003 ± 0.0001 


37.5 


445.9 


0.795 


281.2 


0.0003 ± 0.0001 


42.5 


446.1 


0.795 


281.0 


0.0002 ± 0.0001 


47.5 


448.1 


0.790 


281.8 


0.0002 ± 0.0001 


52.5 


450.5 


0.786 


282.6 


0.0003 ± 0.0001 


57.5 


452.8 


0.781 


283.6 


0.0002 ± 0.0001 



TABLE XVI: Measured cross-section data for QE proton knockout from 16 for 
E beam = 1.643 GcV, 9 pq = 0.0° , and E misB > 25 McV. The p miss bins were 5 
McV/c wide. Cuts were applied to remove the radiative tail from 1 H(e, ep). There 
is an additional 5.8% systematic uncertainty associated with these results. 



Emiss 


< Ul > 


< Q 2 > 


< Pmiss > 


d b a I dujdEpd£l e d$2 p 


(McV) 


(McV) 


(GcV/c) 2 


(McV/c) 


(nb/McV 2 /sr 2 ) 


27.5 


450.0 


0.795 


62.5 


0.0195 ± 0.0013 


32.5 


450.0 


0.795 


62.5 


0.0293 ± 0.0012 


37.5 


450.0 


0.795 


62.5 


0.0483 ± 0.0013 


42.5 


450.0 


0.795 


62.5 


0.0534 ± 0.0014 


47.5 


450.0 


0.795 


62.5 


0.0445 ± 0.0014 


52.5 


450.0 


0.795 


65.0 


0.0263 ± 0.0015 


57.5 


450.0 


0.795 


67.5 


0.0132 ± 0.0021 



32 



TABLE XVII: Measured cross-section data for QE proton knockout from 16 for 
-Ebeam = 1-643 GcV, \0 pq \ = 8.0°, and E miss > 25 McV. The p misg bins were 5 
McV/c wide. There is an additional 6.0% systematic uncertainty associated with 
these results. 




















+8.0 U 






o pq 




-8.0° 




< Ul > 


< 


Q 2 > 


< Pmiss > 


d° a 1 ' dudEpdn^dilp 


< Pmi 


B > 


d a a / ' dLddE p dn e dQ. p 


(McV) 


(MeV) 


(GoV/c) 2 


(McV/c) 


(nb/MeV 2 /sr 2 ) 


(McV/c) 


(nb/MoV 2 /sr 2 ) 


25.5 


436.2 





804 


141.8 


0.0361 


i 





0036 


153 


2 


0.0163 


± 


0.0020 


26.5 


438.0 





802 


142.7 


0.0231 


| 





0031 


152 


2 


0.0130 


± 


0.0018 


27.5 


439.0 





802 


144.1 


0.0193 







0029 


151 


7 


0.0079 


± 


0.0016 


28.5 


439.5 





802 


144.9 


0.0138 


| 





0026 


150 


5 


0.0060 


± 


0.0015 


29.5 


441.2 





800 


146.6 


0.0167 


± 





0026 


149 


7 


0.0117 


± 


0.0016 


30.5 


441.7 


() 


800 


146.8 


0.0209 







0027 


148 


8 


0.0132 


± 


0.0016 


31.5 


442.4 





801 


147.5 


0.0157 


i 





0025 


147 


4 


0.0101 


± 


0.0015 


32.5 


443.6 





799 


148.8 


0.0201 


± 





0026 


147 


2 


0.0075 


± 


0.0014 


33.5 


445.1 





798 


149.9 


0.0185 


| 





0026 


146 


1 


0.0098 


± 


0.0014 


34.5 


445.8 





796 


151.3 


0.0171 


| 





0025 


144 


(i 


0.0138 


± 


0.0015 


35.5 


446.9 





796 


151.1 


0.0178 


| 





0025 


144 


3 


0.0109 


± 


0.0015 


36.5 


447.5 





796 


151.9 


0.0226 


| 





0026 


144 





0.0135 


± 


0.0015 


37.5 


448.9 





795 


153.8 


0.0218 


| 





0025 


142 


7 


0.0128 


± 


0.0015 


38.5 


450.0 





796 


154.5 


0.0245 


| 





0027 


141 


7 


0.0139 


± 


0.0016 


39.5 


450.6 





795 


155.1 


0.0214 


| 





0026 


140 





0.0151 


± 


0.0016 


40.5 


451.5 





795 


156.0 


0.0242 


± 





0027 


141 


1 


0.0128 


± 


0.0015 


41.5 


452.0 





793 


156.6 


0.0264 


[ 





0027 


139 


6 


0.0144 


± 


0.0016 


42.5 


453.7 


() 


792 


158.1 


0.0171 


± 





0025 


137 


() 


0.0173 


± 


0.0017 


43.5 


454.9 





791 


159.3 


0.0175 


± 





0024 


137 


9 


0.0144 


± 


0.0016 


44.5 


456.3 





791 


160.9 


0.0202 







0025 


136 


6 


0.0142 


± 


0.0016 


45.5 


456.6 





791 


160.4 


0.0203 


± 





0025 


135 


5 


0.0166 


± 


0.0016 


46.5 


457.4 





789 


161.3 


0.0136 


| 





0022 


136 


3 


0.0149 


± 


0.0016 


47.5 


459.0 





788 


162.5 


0.0146 


± 





0022 


134 


5 


0.0162 


± 


0.0016 


48.5 


459.9 





788 


163.7 


0.0170 


| 





0023 


133 


7 


0.0127 


± 


0.0015 


49.5 


461.1 





787 


164.6 


0.0144 


| 





0022 


132 


8 


0.0125 


± 


0.0015 


50.5 


462.3 





787 


165.6 


0.0157 


| 


() 


0023 


132 


5 


0.0138 


± 


0.0016 


51.5 


462.8 





786 


166.4 


0.0142 


| 





0021 


131 


9 


0.0110 


± 


0.0015 


52.5 


464.5 





785 


168.0 


0.0079 


± 





0019 


130 


2 


0.0130 


± 


0.0015 


53.5 


464.6 





785 


168.4 


0.0078 


i 





0019 


130 


8 


0.0093 


± 


0.0014 


54.5 


465.5 





784 


169.3 


0.0073 


i 





0018 


130 


1 


0.0112 


± 


0.0014 


55.5 


466.4 





785 


170.1 


0.0096 


i 





0019 


129 


7 


0.0109 


± 


0.0015 


56.5 


467.0 





784 


170.5 


0.0094 


I 





0019 


128 


:j 


0.0069 


± 


0.0013 


57.5 


467.3 





783 


171.1 


0.0066 


± 





0018 


128 


4 


0.0078 


± 


0.0013 


58.5 


468.2 





783 


171.0 


0.0062 


i 





0018 


129 


5 


0.0070 


± 


0.0013 


59.5 


468.3 





782 


171.8 


0.0041 







0016 


127 


2 


0.0097 


± 


0.0014 


60.5 


469.2 





782 


173.1 


0.0100 


± 





0019 


127 


8 


0.0070 


± 


0.0013 


61.5 


469.3 





782 


172.9 


0.0060 







0018 


127 


1 


0.0036 


± 


0.0012 


62.5 


470.6 


() 


782 


173.9 


0.0051 


± 





0018 


127 


3 


0.0085 


± 


0.0014 


63.5 


469.9 





782 


173.6 


0.0062 


i 





0017 


126 


8 


0.0040 


± 


0.0012 


64.5 


471.0 





781 


174.7 


0.0071 


i 





0019 


126 


9 


0.0065 


± 


0.0013 


65.5 


471.2 





780 


174.5 


0.0045 


± 





0018 


125 


8 


0.0043 


± 


0.0013 


66.5 


472.1 





779 


175.5 


0.0077 


| 





0019 


125 


9 


0.0056 


± 


0.0013 


67.5 


472.6 





780 


175.9 


0.0049 


± 





0019 


125 


5 


0.0019 


± 


0.0012 


68.5 


473.4 





779 


177.2 


0.0061 


| 


() 


0019 


126 


1 


0.0049 


± 


0.0013 


69.5 


473.3 





779 


177.0 


0.0071 


| 





0019 


124 


() 


0.0022 


± 


0.0012 


70.5 


474.4 





778 


177.9 


0.0020 


| 


() 


0017 


125 


2 


0.0056 


± 


0.0013 


71.5 


474.5 





779 


178.7 


0.0033 







0017 


124 


8 


0.0034 


± 


0.0013 


72.5 


474.9 





778 


178.4 


0.0052 


± 





0018 


125 


2 


0.0066 


± 


0.0015 


73.5 


475.1 





778 


179.0 


0.0065 


b 





0020 


125 


1 


0.0046 


± 


0.0014 


74.5 


476.3 





777 


179.6 


0.0027 


i 





0018 


124 


4 


0.0032 


± 


0.0013 


75.5 


476.2 





776 


180.3 


0.0026 


1 





0017 


124 


4 


0.0046 


± 


0.0014 


76.5 


477.1 





776 


181.4 


0.0061 


1 





0021 


122 


3 


0.0005 


± 


0.0012 


77.5 


477.6 





776 


181.7 


0.0034 


b 





0019 


124 


1 


0.0016 


± 


0.0013 


78.5 


478.2 





777 


182.7 


0.0056 


1 





0021 


123 


8 


0.0013 


± 


0.0014 


79.5 


478.9 





775 


183.1 


0.0056 


■1 





0021 


123 


4 


0.0054 


± 


0.0016 


80.5 


479.1 





774 


183.3 


0.0059 


b 





0023 


122 





0.0040 


± 


0.0016 


81.5 


479.6 





776 


184.1 


0.0008 


t 





0018 


122 


7 


0.0039 


± 


0.0017 


82.5 


480.1 





775 


185.1 


0.0007 


± 





0018 


122 


8 


0.0029 


± 


0.0016 


83.5 


480.6 





776 


185.2 


0.0039 


t 





0021 


121 


8 


0.0022 


± 


0.0016 


84.5 


481.1 





773 


185.8 


0.0058 


1 





0025 


122 


8 


0.0018 


± 


0.0016 


85.5 


481.6 





775 


186.3 


0.0028 


± 





0023 


122 


2 


0.0039 


± 


0.0019 


86.5 


482.1 





775 


187.1 


0.0065 


1 





0027 


122 


9 


0.0003 


± 


0.0018 


87.5 


482.8 





772 


187.6 


0.0049 


1 





0030 


123 


1 


0.0029 


± 


0.0021 
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TABLE XVIII: Measured cross-section data for QE proton knockout from 16 O for 
Ebeam = 2.442 GeV, 9 pq = 0.0° , and E miss > 25 MeV. The p miss bins were 5 McV/c 
wide. Cuts were applied to remove the radiative tail from 1 H(e, ep). There is an 
additional 5.8% systematic uncertainty associated with these results. 





< ui > 


< Q 2 > 


< Pmiss > 




(MeV) 


(MeV) 


(GoV/c) 2 


(McV/c) 


(nb/McV 2 /sr 2 ) 


27.5 


450.0 


0.795 


60.0 


0.0552 ± 0.0042 


32.5 


450.0 


0.795 


60.0 


0.0890 ± 0.0045 


37.5 


450.0 


0.795 


60.0 


0.1387 ± 0.0050 


42.5 


450.0 


0.795 


60.0 


0.1580 ± 0.0058 


47.5 


450.0 


0.795 


62.5 


0.1348 ± 0.0062 


52.5 


450.0 


0.795 


65.0 


0.0756 ± 0.0063 


57.5 


450.0 


0.795 


67.5 


0.0402 ± 0.0084 



TABLE XIX: Measured cross-section data for QE proton knockout from ls O for 
E beam = 2.442 GeV, 9 pq = -2.5°, and E miss > 25 MeV. The p misB bins were 5 
McV/c wide. There is an additional 5.9% systematic uncertainty associated with 
these results. 



Bmiss 


<u > 


< 


Q 2 > 


< Pmiss > 


d° a I ' dujdE p dQ. e d{1 v 


(MeV) 


(MeV) 


(GcV/c) 2 


(McV/c) 


(nb/McV 2 /sr 2 ) 


25.5 


445.0 





823 


46.5 





0330 


± 


0.0182 


26.5 


446.5 





822 


46.8 





0661 


± 


0.0185 


27.5 


449.5 





821 


46.7 





0388 


± 


0.0173 


28.5 


449.3 





820 


45.3 





0590 


± 


0.0176 


29.5 


450.4 





823 


46.2 





0955 


± 


0.0192 


30.5 


450.8 





823 


45.0 





0589 


± 


0.0181 


31.5 


451.9 





821 


45.4 





0733 


± 


0.0180 


32.5 


452.1 





823 


45.0 





0810 


± 


0.0179 


33.5 


453.3 





822 


45.6 





0944 


± 


0.0192 


34.5 


454.5 





821 


45.2 





1375 


± 


0.0211 


35.5 


455.1 





821 


46.1 





0984 


± 


0.0209 


36.5 


454.9 





822 


46.1 





1950 


± 


0.0227 


37.5 


457.3 





821 


45.6 





1566 


± 


0.0232 


38.5 


458.4 





820 


45.8 





1937 


± 


0.0232 


39.5 


458.4 





820 


46.1 





2052 


± 


0.0251 


40.5 


459.3 





821 


46.5 





1841 


± 


0.0236 


41.5 


459.9 





819 


45.8 





2354 


± 


0.0252 


42.5 


461.3 





818 


45.9 





2401 


± 


0.0260 


43.5 


462.7 





817 


45.6 





1791 


± 


0.0247 


44.5 


462.6 





818 


46.6 





2049 


± 


0.0245 


45.5 


464.0 





816 


45.9 





1401 


± 


0.0223 


46.5 


464.7 





816 


46.9 





1721 


± 


0.0224 


47.5 


467.9 





817 


45.7 





1457 


± 


0.0217 


48.5 


467.7 





816 


45.9 





1323 


± 


0.0211 


49.5 


468.2 





815 


45.8 





1469 


± 


0.0211 


50.5 


471.1 





814 


45.3 





0893 


± 


0.0189 


51.5 


471.0 





814 


45.9 





1008 


± 


0.0185 


52.5 


473.1 





814 


46.3 





0866 


± 


0.0175 


53.5 


472.1 





813 


45.8 





1037 


± 


0.0183 


54.5 


474.6 





813 


45.7 





0607 


± 


0.0165 


55.5 


474.5 





812 


45.4 





0766 


± 


0.0163 


56.5 


473.8 





810 


46.8 





0617 


± 


0.0155 


57.5 


476.5 





810 


45.8 





0606 


± 


0.0155 


58.5 


477.4 





811 


46.6 





0717 


± 


0.0158 


59.5 


479.4 





809 


45.1 





0837 


± 


0.0168 


60.5 


481.5 





810 


45.5 





0298 


± 


0.0140 


61.5 


480.8 





807 


46.2 





0748 


± 


0.0156 


62.5 


480.6 





807 


47.7 





0165 


± 


0.0125 


63.5 


480.7 





805 


47.8 





0240 


± 


0.0118 


64.5 


481.1 





805 


47.6 





0135 


± 


0.0109 


65.5 


484.1 





805 


46.5 





0471 


± 


0.0127 


66.5 


482.4 





804 


46.6 





0195 


± 


0.0118 


67.5 


485.9 





804 


47.4 





0297 


± 


0.0121 


68.5 


484.4 





803 


48.7 





0491 


± 


0.0134 


69.5 


490.6 





802 


46.5 





0118 


± 


0.0116 


70.5 


489.5 





800 


47.6 





0021 


± 


0.0103 
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TABLE XX: Measured cross-section data for QE proton knockout from ie O for 
Ebeam = 2.442 GcV, \9 pq \ = 8.0°, and E miss > 25 McV. The p miss bins were 5 
McV/c wide. There is an additional 6.0% systematic uncertainty associated with 
these results. 
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TABLE XXI: Measured cross-section data for QE proton knockout from 16 O for 
Ebeam = 2.442 GcV, \6 pq \ = 16.0°, and E miss > 25 McV. The p miss bins were 5 
McV/c wide. There is an additional 5.8% systematic uncertainty associated with 
these results. 
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TABLE XXII: Measured cross-section data for QE proton knockout from 16 for 
-Ebcam = 2.442 GeV, \0 pq \ = 20.0°, and E miss > 25 MeV. The p miss bins were 5 
McV/c wide. There is an additional 5.9% systematic uncertainty associated with 
these results. 
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2. Effective response functions and asymmetries 



a. lp-shell 



Effective response functions and asymmetries for QE proton knockout from 16 arc presented in Tables XXIII. 
XXIV, and XXV. 
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TABLE XXIII: Effective response functions Rl and Rt for QE proton knockout 
from the lp-shell of ie O for < T p > = 427 McV. The p miss bins were 20 McV/c wide. 
Cuts were applied to remove the radiative tail from 1 H(e, ep) such that < p m i ss > — 
52.5 MeV/c in each case. Note that the u> acceptance of the three measurements was 
shifted by about 3% in order to keep the phase-space acceptance flat. The resulting 
change in the domain was taken as an additional 4% systematic uncertainty. 



lpi^-statc 


lp 3/2 -statc 


< Pmiss > 

(McV/c) 


Rl (fm 3 ) 
(stat) (sys) 


Rt (fm 3 ) 
(stat) (sys) 


< Pmiss > 

(MeV/c) 


Rl (fm 3 ) 
(stat) (sys) 


Rt (fm 3 ) 
(stat) (sys) 


52.5 


1.82 ± 1.17 ± 0.59 


7.58 ± 1.42 ± 0.72 


52.5 


2.35 ± 1.18 ± 0.73 


9.40 ± 1.39 ± 0.85 



TABLE XXIV: Effective response functions Rl+tt and Rt for QE proton knockout 
from the lp-shell of le O for < Q 2 > = 0.800 (GeV/c) 2 , < u > = 436 McV, and 
< T p > = 427 McV. The p miee bins were 20 MeV/c wide. 



lpl/2 -ST;a te 


lp3/2-statc 


< Pmiss > 

(McV/c) 


Rl+tt (fm 3 ) 
(stat) (sys) 


Rt (fm 3 ) 
(stat) (sys) 


< Pmiss > 

(MeV/c) 


Rl + tt (fm 3 ) 
(stat) (sys) 


Rt (fm 3 ) 
(stat) (sys) 


149.0 
279.0 


0.56 ± 0.49 ± 0.12 
0.01 ± 0.03 ± 0.01 


6.08 ± 0.61 ± 0.24 
0.12 ± 0.04 ± 0.01 


149.0 
276.0 


2.20 ± 0.75 ± 0.30 
0.20 ± 0.06 ± 0.03 


10.35 ± 1.04 ± 0.43 
0.29 ± 0.08 ± 0.01 



TABLE XXV: The asymmetry Alt and effective response functions Rlt for QE 
proton knockout from the lp-shell of le O for < Q 2 > = 0.800 (GeV/c) 2 , < ui > = 
436 McV, and < T p > = 427 MeV. The p misB bins were 20 MeV/c wide. Save for 
the data labelled [*] which were obtained at Ebeam — 1.643 GeV, the beam energy 
was 2.442 GeV. 



lpi/2 _s tatc 


lP3/2 _s tatc 


< Pmiss > 

(MeV/c) 


Alt 
(stat) (sys) 


Rlt (fm 3 ) 
(stat) (sys) 


< Pmiss > 

(MeV/c) 


Alt 
(stat) (sys) 


Rlt (fm 3 ) 
(stat) (sys) 


60.0 
148.0 [*] 
149.0 
279.0 
345.0 


0.02 ± 0.02 ± 0.02 
-0.25 ± 0.02 ± 0.03 
-0.23 ± 0.02 ± 0.03 
-0.36 ± 0.08 ± 0.04 
-0.13 ± 0.22 ± 0.05 


0.117 ± 0.134 ± 0.037 
-1.198 ± 0.235 ± 0.085 
-0.999 ± 0.066 ± 0.077 
-0.029 ± 0.007 ± 0.002 
-0.002 ± 0.003 ± 0.001 


60.0 
147.0 [*] 
148.0 
276.0 
345.0 


-0.08 ± 0.02 ± 0.02 
-0.31 ± 0.04 ± 0.03 
-0.31 ± 0.01 ± 0.03 
-0.69 ± 0.04 ± 0.04 
-0.39 ± 0.08 ± 0.05 


-0.754 ± 0.165 ± 0.084 
-2.820 ± 0.292 ± 0.183 
-2.560 ± 0.096 ± 0.173 
-0.250 ± 0.013 ± 0.019 
-0.015 ± 0.003 ± 0.001 



b. Higher missing energies 



Effective response functions for QE proton knockout from 16 for E m - lss > 25 MeV are presented in Tables XXVI 
- XXVIII. 

TABLE XXVI: Effective response functions Rl, and Rt for QE proton knockout 
from ie O at 9 pq = 0° for E miss > 25 MeV. The p mias bins were 5 McV/c wide. 



< Bmiss > 


Rl (fm 3 /McV) 


R T (fm 3 /MeV) 


(McV) 


(stat) (sys) 


(stat) (sys) 


27.5 


0.067 ± 0.050 ± 0.020 


0.203 ± 0.065 ± 0.021 


32.5 


0.075 ± 0.044 ± 0.030 


0.346 ± 0.056 ± 0.035 


37.5 


0.141 ± 0.040 ± 0.052 


0.523 ± 0.049 ± 0.051 


42.5 


0.138 ± 0.037 ± 0.063 


0.600 ± 0.044 ± 0.066 


47.5 


0.100 ± 0.039 ± 0.054 


0.517 ± 0.047 ± 0.051 


52.5 


0.054 ± 0.044 ± 0.049 


0.304 ± 0.054 ± 0.028 


57.5 


-0.047 ± 0.111 ± 0.030 


0.253 ± 0.145 ± 0.049 



TABLE XXVII: Effective response functions Rl+tt and Rt for QE proton knock- 
out from le O at 9 pq = 8° and 8 pq = 16° for E m iss > 25 McV. The p m iss bins were 
5 MeV/c wide. 





e pg = 16 


< Bmiss > 

(MeV) 


Rl+tt (fm 3 /McV) 
(stat) (sys) 


Rt (fm 3 /MeV) 
(stat) (sys) 


Rl + tt (fm 3 /MeV) 
(stat) (sys) 


R T (fm 3 /McV) 
(stat) (sys) 


27.5 
32.5 
37.5 
42.5 
47.5 
52.5 
57.5 


0.080 ± 0.023 ± 0.010 
0.072 ± 0.021 ± 0.008 
0.075 ± 0.021 ± 0.008 
0.087 ± 0.021 ± 0.013 
0.050 ± 0.022 ± 0.008 
0.001 ± 0.022 ± 0.006 
-0.002 ± 0.024 ± 0.004 


0.137 ± 0.029 ± 0.008 
0.133 ± 0.026 ± 0.008 
0.162 ± 0.027 ± 0.010 
0.164 ± 0.026 ± 0.010 
0.172 ± 0.028 ± 0.010 
0.137 ± 0.027 ± 0.009 
0.010 ± 0.030 ± 0.006 


0.002 ± 0.004 ± 0.003 
-0.001 ± 0.004 ± 0.001 
0.001 ± 0.004 ± 0.001 
0.005 ± 0.004 ± 0.001 
0.001 ± 0.004 ± 0.001 
0.005 ± 0.005 ± 0.001 
0.008 ± 0.007 ± 0.001 


0.027 ± 0.005 ± 0.005 
0.026 ± 0.005 ± 0.001 
0.018 ± 0.004 ± 0.001 
0.014 ± 0.004 ± 0.001 
0.017 ± 0.005 ± 0.001 
0.018 ± 0.006 ± 0.002 
0.008 ± 0.007 ± 0.001 
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TABLE XXVIII: Effective response functions Rlt for QE proton knockout from 
ie O at pq = 8° and 6 vq = 16° for E miss > 25 MeV. The p miss bins were 5 McV/c 
wide. Note that the data presented for 6 vq — 8° represent the average of the data 
obtained at E hcaTn = 1.643 GcV and -B beam = 2.442 GeV. 



Opq— 8 


e pq = 16" 


< -E miB3 > 
(MoV) 


R LT (fm'VMoV) 
(stat) (sys) 


R LT (fm J /McV) 
(stat) (sys) 


27.5 


-0.070 ± 0.006 ± 0.005 


-0.009 ± 0.001 ± 0.001 


32.5 


-0.050 ± 0.005 ± 0.004 


-0.007 ± 0.001 ± 0.001 


37.5 


-0.076 ± 0.005 ± 0.010 


-0.007 ± 0.001 ± 0.002 


42.5 


-0.083 ± 0.006 ± 0.039 


-0.007 ± 0.001 ± 0.001 


47.5 


-0.077 ± 0.006 ± 0.010 


-0.004 ± 0.002 ± 0.001 


52.5 


-0.036 ± 0.006 ± 0.015 


-0.004 ± 0.002 ± 0.001 


57.5 


-0.036 ± 0.006 ± 0.005 


-0.004 ± 0.003 ± 0.001 



APPENDIX B: A 'DIP'-REGION 
INVESTIGATION 

A small portion of the beam time allocated to the mea- 
surement discussed in the main body of this article was 
used for an exploratory investigation of the 'dip' located 
in the energy-transfer region between the QE peak and 
the A(1232)-resonance. For this investigation, -Ebcam = 
1.643 GeV was employed, and the HRS e position and cen- 
tral momentum were fixed at 9 e = 37.17° and p e = 1056 
MeV/c, respectively. This resulted in q ps 1.026 GeV/c, ui 
ps 589 MeV, and Q 2 ps 0.706 (GeV/c) 2 [148]. The HRS h 
was then positioned at Oh = 38.45° (9 pq = 0°) and its cen- 
tral momentum varied from 828 McV/c to 1190 MeV/c 
in five steps of Ap p ps 70 McV/c per step. These momen- 
tum settings were close enough to each other that there 
was adequate acceptance overlap between them to allow 
for radiative corrections to be performed. The configura- 
tion of the experimental apparatus and data-acquisition 
system was identical in all aspects to that used for the 
QE measurement. The data analysis was also identical 
to that performed on the QE data, save for an additional 
cut to remove H(e, e'p)ir° events. 

Fig. 25 shows the measured cross-section data for the 
dip region as a function of i? m i S s compared to calculations 
by the Ghent Group for -Ebcam = 1-643 GeV (see Table 
XXIX) . The dashed curve is the bare romea calculation 
for proton knockout from the lsi/2-state of 16 and the 
solid curve is the same calculation including the effects of 
MEC and IC (see the main text of this article for further 
details). A normalization factor of 1.0 was employed for 
these calculations. The dashed-dotted curve illustrates 
the calculated (e, e'pN) contribution. 

In contrast to the QE energy region, the bare calcu- 
lation actually overestimated the ls^-state strength in 
these kinematics. Also in contrast to the QE energy re- 



gion, the inclusion of MEC and IC decrease the mag- 
nitude of the calculated cross section and improve the 
agreement. Finally, while the (e,e'pN) calculations have 
the measured flat shape for E m i ss > 100 MeV, they are 
twice as large as the cross-section data. 
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FIG. 25: Data from this work together with calculations by 
the Ghent Group for the .E m i ss -dependence of the cross-section 
data obtained in dip-region kinematics for -Ebcam = 1.643 
GeV. Uncertainties are statistical and, on average, there is an 
additional ±5.9% systematic uncertainty (see Table XXIX) 
associated with the data. 
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TABLE XXIX: Measured cross-section data for proton knock- 
out from 16 O for -Eb cam = 1.643 GeV in dip- region kinematics 
for -Emias > 25 MeV. There is an additional 5.9% systematic 
uncertainty associated with these results. 



B mia , 


< Pmiss > 


d G a j dujdE p dQ. e di~l v 


(McV) 


(MoV/c) 


(nb/M G V 2 /sr 2 ) 


27.5 


163.6 


0.0115 


± 


0.0012 


32.5 


152.1 


0.0158 


± 


0.0012 


37.5 


138.8 


0.0235 


± 


0.0013 


42.5 


132.8 


0.0316 


± 


0.0015 


47.5 


128.3 


0.0262 


± 


0.0013 


52.5 


119.4 


0.0184 


± 


0.0011 


57.5 


112.1 


0.0160 


± 


0.0011 


62.5 


107.5 


0.0107 


± 


0.0009 


67.5 


95.0 


0.0078 


± 


0.0009 


72.5 


95.0 


0.0079 


± 


0.0008 


77.5 


95.0 


0.0062 


± 


0.0007 


82.5 


90.0 


0.0045 


± 


0.0006 


87.5 


84.0 


0.0028 


± 


0.0005 


92.5 


76.0 


0.0025 


± 


0.0005 


97.5 


71.7 


0.0030 


± 


0.0005 


102.5 


68.0 


0.0025 


± 


0.0005 


107.5 


60.0 


0.0015 


± 


0.0004 


112.5 


61.4 


0.0022 


± 


0.0004 


117.5 


58.0 


0.0013 


± 


0.0004 


122.5 


51.0 


0.0006 


± 


0.0004 


127.5 


47.5 


0.0010 


± 


0.0004 


132.5 


50.0 


0.0015 


± 


0.0005 


137.5 


60.0 


0.0007 


± 


0.0005 


142.5 


45.0 


0.0012 


± 


0.0006 


147.5 


39.0 


0.0002 


± 


0.0004 


152.5 


39.3 


0.0014 


± 


0.0004 


157.5 


48.8 


0.0012 


± 


0.0004 


162.5 


51.7 


0.0009 


± 


0.0004 


167.5 


55.0 


0.0019 


± 


0.0004 


172.5 


62.0 


0.0016 


± 


0.0004 


177.5 


65.9 


0.0015 


± 


0.0004 


177.5 


65.9 


0.0015 


± 


0.0004 


182.5 


69.0 


0.0028 


± 


0.0005 


187.5 


77.2 


0.0023 


± 


0.0005 


192.5 


94.0 


0.0017 


± 


0.0004 


197.5 


97.7 


0.0016 


± 


0.0005 


202.5 


97.7 


0.0020 


± 


0.0005 


207.5 


103.0 


0.0015 


± 


0.0005 


212.5 


114.1 


0.0014 


± 


0.0005 


217.5 


118.3 


0.0009 


± 


0.0004 


222.5 


122.3 


0019 




o nnn5 


227.5 


125.9 


0.0020 


± 


0.0005 


232.5 


135.0 


0.0013 


± 


0.0005 


237.5 


150.5 


0.0014 


± 


0.0006 


242.5 


157.7 


0.0011 


± 


0.0005 


247.5 


157.7 


0.0004 


± 


0.0004 


252.5 


167.0 


0.0018 


± 


0.0005 


257.5 


170.0 


0.0013 


± 


0.0005 


262.5 


173.2 


0.0018 


± 


0.0006 


267.5 


179.5 


0.0022 


± 


0.0006 


272.5 


185.9 


0.0008 


± 


0.0005 


277.5 


188.0 


0.0012 


± 


0.0006 


282.5 


210.0 


0.0006 


± 


0.0005 


287.5 


200.0 


0.0017 


± 


0.0006 


292.5 


200.0 


0.0004 


± 


0.0006 


297.5 


205.0 


0.0009 


± 


0.0006 
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[133] Kinematically, an electron scattered through angle 8 e 
transfers momentum q and energy lu with Q 2 = q 2 — lu 2 . 
The ejected proton has mass m p , momentum p p , en- 
ergy E p , and kinetic energy T p . In QE kinematics, 
lu % Q 2 /2m p . The cross section is typically measured 
as a function of missing energy -E m i ss = lu — T p — Tb 
and missing momentum p m i sa = \q — p p \- Tb is the ki- 
netic energy of the residual nucleus. The lab polar an- 
gle between the ejected proton and virtual photon is pq 
and the azimuthal angle is cj>- B vq > 0° corresponds to 
<f> — 180°, 8 P > q , and +p m i ss . 6 pq < 0° corresponds to 
4> = 0° , 6 P < 8 q , and — p m i ss . 

[134] In the One-Photon Exchange Approximation, the un- 
polarized (e, e'p) cross section can be expressed as the 
sum of four independent response functions: Rl (longi- 
tudinal), Rt (transverse), Rlt (longitudinal-transverse 
interference, and Rtt (transverse-transverse interfer- 
ence). See also Eq. (4). 



[135] In the nonrelativistic Plane- Wave Impulse Approxima- 
tion, the transverse amplitude in the Rlt response 
is uniquely determined by the convection current. At 
higher Q 2 , it is well-known that the convection current 
yields small matrix elements. As a result, the nonrela- 
tivistic Impulse Approximation (ia) contributions which 
dominate Rl and Rt are suppressed in Rlt (and thus 
Alt)- Hence, these observables are particularly sensi- 
tive to any mechanisms beyond the IA, such as channel 
coupling and relativistic and two-body current mecha- 
nisms [105]. 

[136] Alt = ^ ( ( $_io°)+^(0=i8O°) ■ Alt is a Particularly useful 
quantity for experimentalists because it is systemati- 
cally much less challenging to extract than either an 
absolute cross section or an effective response function. 

[137] ls-shell nucleons are generally knocked out from high- 
density regions of the target nucleus. In these high- 
density regions, the IA is expected to be less valid 
than for knockout from the valence p-shell states lying 
near the surface. In this region of 'less-valid' IA, size- 
able contributions to the s-shell cross-section data arise 
from two-nucleon current contributions stemming from 
MEC and IC. In addition to affecting the single-nucleon 
knockout cross section, the two-nucleon currents can re- 
sult in substantial multi-nucleon knockout contributions 
to the higher _B m i ss continuum cross section [105]. 

[138] The transverse- longitudinal difference is St — Sl, where 
Sx = o~MottVx Rx/o-^p, and X e {T,L}. a* p represents 
components of the off-shell ep cross section and may be 
calculated using the CCl, CC2, or CC3 prescriptions of 
de Forest [83]. 

[139] When necessary, the differential dependencies of the 
measured cross-section data were changed to match 
those employed in the theoretical calculations. The pris- 
tine detection volume AVb(£ , m i S s, Pmiss, lu, Q 2 ) was 
changed to a weighted detection volume by weighting 
each of the trials with the appropriate Jacobian(s). 

[140] The difference between cross-section data averaged over 
the reduced spectrometer acceptances and calculated for 
a small region of the central kinematics was no more 
than 1%. Thus, the finite acceptance of the spectrome- 
ters was not an issue. 

[141] This Jacobian is given by = §^ + ^|f-, where 

Eb = y/p% + m%. 
[142] The phase-space factor K is given in Eq. (6) , while CMott 

= T 'i C ° s d[ 6 n ° ^1 ■ The dimensionless kinematic factors are 

4sf sm 4 (0 e /2) 

as follows: vl = yr, vt — -^p- + tan 2 (# e /2), vlt = 

f^f^+tan 2 (0 e /2), and v TT = _fr- 

[143] The accuracy of the effective response-function sepa- 
ration depends on precisely matching the values of q 
and lu at each of the different kinematic settings. This 
precise matching was achieved by measuring 1 H(e,ep) 
with a pinhole collimator (in practice, the central hole 
of the sieve-slit collimator) placed in front of the HRS e . 
The proton momentum was thus q. The 1 H(e,ep) pro- 
ton momentum peak was determined to Ap/p — 1.5 x 
10 -4 , which allowed for an identical matching of Aq/q 
between the different kinematic settings. 

[144] Rl+tt = Rl + ^Rtt- 

[145] Note that the calculations in Ref. [2] using LEA neglected 
the (S — V) term and replaced the momentum in Eq. 
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(29b) by its asymptotic value, an approach later called 
EMA-noSV, where EMA denotes the Effective Momen- 
tum Approximation. 

[146] Strictly speaking, the effective longitudinal response 
function Rl could not be separated from the quasiper- 
pendicular kinematics data. However, since both Kelly 
and Udfas et al. calculate the term ^^-Rtt to be < 
10% of Rl+tt in these kinematics, Rl and Rl+tt re- 
sponses are both presented on the same plot. 

[147] Difficulties associated with the reaction mechanism rel- 
atively far from QE kinematics may also be partly re- 



ponsible for the anomalously low normalization factors 
for these data. 

[148] The quantity y (which is the minimum value of the 
initial momentum of the nucleon) is generally used 
to label non-QE kinematics. According to Day et al. 
[132], y = {(m A + uj) ^/ h 2 - m 2 B W 2 - qA]/W 2 , where 
W = ,/(mA+a/) 2 -q 2 and A = (m| — m 2 N + W 2 )/2. 
In QE kinematics, y — 0. y — 0.16 in these dip-region 
kinematics. 



